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ABSTRACT 


A  mathematical  model  is  developed  to  analyze  the  seismic 
response  of  a  drydocked  vessel  in  three  degrees  of  freedom; 
vessel  rotation  about  its  keel  and  vessel  horizontal  and 
vertical  translations  relative  to  the  drydock  cradle.   Data 
from  eleven  actual  vessel-drydock  systems  and  the  time 
acceleration  history  of  an  earthquake  are  implemented  to 
predict  vessel  three  degree  of  freedom  response  during  an 
earthquake.   Vessel  seismic  response  and  the  resultant  drydock 
forces  are  compared  to  vessel-drydock  system  failure  criteria 
to  determine  stability  of  the  system  during  earthquakes.   The 
three  degree  of  freedom  vessel  response  model  is  compared  to  a 
one  and  two  degree  of 'freedom  vessel  response  models  and  a 
model  in  which  seismic  loading  is  simulated  by  a  single  static 
force.   The  three  degree  of  freedom  vessel  motion  is  shown  to 
be  the  most  accurate  method  for  analyzing  vessel-drydock 
system  failure  criteria. 
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1 . 0   INTRODUCTION 

1.1.   Background 

In  the  design  of  a  drydock  for  military  vessels,  an 
important  factor  is  the  environmental-generated  forces  that 
are  found  in  the  drydock  cradle  which  supports  the  ship. 
Earthquakes  and  the  wind  are  considered  the  two  most 
significant  environmental  loading  factors.   This  thesis 
explores  the  seismic  loading  aspect  of  the  problem. 

Regardless  of  its  location,  any  U.S.  Navy  ship  that  is 
nuclear  powered  and  in  drydock  is  considered  to  be  in  a  high 
seismic  risk  area.   Currently  the  Navy  defines  a  high  seismic 
risk  area  as  one  in  which  earthquake  forces  be  approximated  by 
a  steady  horizontal  force  of  0.2  g  times  the  vessel  mass 
acting  at  the  center  of  gravity.   The  drydock  can  be  analyzed 
to  see  if  it  can  withstand  this  quasi-static  loading. 

1.2   Previous  Work 

In  1981  B.V.  Viscomi  studied  the  seismic  response  of  a 
drydocked  submarine  using  a  single  degree  of  freedom  model 
[14].   This  analysis  assumed  that  the  vessel  was  allowed  to 
rotate  about  the  keel.   Viscomi  used  the  time  history 
acceleration  record  of  the  North-South  component  of  the  194  6 
El  Centro,  CA  Earthquake,  California  Institute  of  Technology 
processing  scheme.   The  vessel  was  analyzed  for  system  failure 
(i.e.  the  vessel  lifting  off  of  a  row  of  side  blocks). 

A  thesis  by  C.F.  Barker  in  1985  [1]  used  a  two  degree  of 
freedom  model  to  study  seismic  response.   Employing  the  1946 
El  Centro  Earthquake,  Massachusetts  Institute  of  Technology 
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standard  processing  scheme,  Barker  analyzed  rotation  about  the 
keel  plus  horizontal  relative  translation  between  the  vessel 
and  the  drydock  cradle  perpendicular  to  the  vessel  longitudal 
axis.   The  two  degree  of  freedom  model  permitted  the  analysis 
of  drydock  cradle  failure:   block  sliding  and  tipping,  which 
were  not  possible  to  analyze  in  Viscomi's  one  degree  of 
freedom  model,  along  with  block  liftoff. 

1.3  Contributions  of  This  Thesis 

This  thesis  develops  a  three  degree  of  freedom  model  of  a 
vessel  in  drydock.   Motion  will  be  permitted  in  three 
directions:   rotation  about  the  keel  and  relative 
translations,  horizontal  and  vertical,  between  the  vessel  and 
the  drydock  cradle.   The  three  degree  of  freedom  model  will 
make  possible  the  analysis  of  drydock  cradle  failure: 
additional  block  crushing  forces  due  to  vertical  motion,  which 
was  not  conceivable  in  the  previous  one  and  two  degree  of 
freedom  models  developed  by  Viscomi  and  Barker  respectively, 
along  with  block  sliding,  tipping  and  liftoff. 

The  seismic  input  for  the  three  degree  of  freedom  model 
is  the  time  acceleration  history  of  the  North-South  component 
of  the  1946  El  Centro,  CA  Earthquake,  Massachusetts  Institute 
of  Technology  standard  processing  scheme.   The  MIT  version  of 
the  El  Centro  Earthquake  acceleration  record  is  described  in 
reference  [9]  ,  and  is  displayed  graphically  in  Figure  1.1. 

1.4  Outline  of  This  Thesis 

Section  2  of  this  thesis  describes  the  vessel-drydock 
system  for  the  reader  and  introduces  terms  which  are  used 
throughout  the  thesis.   Section  3  describes  the  various 
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failure  mechanisms  found  in  the  vessel-drydock  system  which 
can  occur. 

Section  4  examines  the  modelling  of  the  vessel-drydock 
system  under  seismic  loading.   Section   4.1  describes  the 
quasi-static  loading  method,  and  Section  4.2  develops  the  one 
degree  of  freedom  equation  of  vessel  motion  model.   Section 
4 . 3  develops  the  two  degree  of  freedom  equations  of  vessel 
motion  model.   Finally,  Section  4.4  develops  the  three  degree 
of  freedom  equations  of  motion  model. 

Section  5  discusses  the  system  parameters  that  are 
required  for  the  mathematical  model  found  in  Section  4 . 
Included  in  this  Section  is  the  modelling  of  block  stiffness. 
Eleven  typical  vessel-drydock  system  configurations  are 
studied  to  be  implemented  into  Section  6. 

Section  6  examines  methods  of  evaluating  the  one,  two  and 
three  degree  of  freedom  models  found  in  Section  4.   Section 
6.1  evaluates  the  linear  equations  of  motion  models  using 
modal  analysis  method.   Section  6.2  evaluates  both  linear  and 
non-linear  equations  of  motion  models  using  a  fourth-order 
Runge-Kutta  numerical  analysis  scheme.  Section  6.3  explains 
the  response  spectrum  method  of  determining  the  maximum  value 
of  vessel  seismic  response.   Section  6.4  determines  system 
response  using  quasi-static  force  method. 

Section  7  describes  the  development  and  testing  of  the 
one,  two  and  three  degree  of  freedom  vessel  motion  computer 
programs,  and  contains  the  results  for  a  seismic  analysis  of 
several  actual  vessel-drydock  systems. 

Section  8  summarizes  the  response  predictions  for  several 
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actual  vessel -drydock  systems  obtained  with  quasi-static,  one 
degree  of  freedom,  two  degree  of  freedom  and  three  degree  of 
freedom  analysis  methods,  and  conclusions  are  drawn. 

Appendix  1  is  an  example  listing  of  the  fourth-order 
Runge-Kutta  computer  program.   Appendix  2  is  an  example 
listing  of  the  modal  analysis  computer  program.   Finally, 
Appendix  3  describes  the  modal  analysis  of  the  two  and  three 
degree  of  freedom  vessel-drydock  systems  and  predicts  maximum 
system  response  using  the  response  spectrum  method  of  analysis 
and  modal  participation  factors. 
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FIGURE  1.1 
Acceleration  Time  History,  El  Centro  Earthquake,  1946 
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2.0   THE  VESSEL-DRYDOCK  SYSTEM 

2.1  General  Description 

In  order  to  perform  inspection,  maintentance  and  repair 
on  the  outer  hull  of  a  ship,  it  is  necessary  to  place  the 
vessel  in  a  drydock.   A  drydock  is  a  concrete  encasement  built 
into  the  earth.   At  one  end  of  the  drydock  is  a  flood  gate  to 
allow  the  entry  and  departure  of  the  vessel.   A  cradle,  formed 
by  drydock  blocks,  is  built  upon  the  drydock  floor  in  order  to 
support  the  ship  once  the  water  is  pumped  out  of  the 
encasement.   The  drydock  must  be  able  to  hold  the  weight  of 
the  vessel  and  cradle. 

Due  to  the  various  hull  configurations  associated  with 
the  different  classes  of  ships,  the  cradle  component  of  the 
drydock-cradle-vessel  system  is  varied  to  suit  a  particular 
docking  situation.   The  cradle  accomplishes  the  transfer  of 
vessel  weight  to  the  drydock  floor,  provides  stable  support  to 
the  vessel  and  allows  access  to  the  vessel  hull.   This  cradle 
is  usually  constructed  of  timber,  concrete,  or  a  timber- 
concrete  composite. 

2 . 2  System  Component 

The  primary  components  of  the  drydock-cradle-vessel 
system  are  shown  in  Figure  2.1.   This  section  defines  terms, 
related  to  docking,  used  in  this  thesis  to  describe  the  system 
components . 
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1.  Blocks  -  The  units,  consisting  of  timber,  concrete,  steel 
and  other  materials,  which  together  make  up  the  cradle. 

2 .  Cradle  -  A  framework  of  blocks  which  supports  the  vessel 
when  the  drydock  is  dewatered. 

3.  ICG  -  The  distance  between  the  vessel  baseline  (ship's 
keel)  to  its  vertical  center  of  gravity. 

4.  Keel  Blocks  -  The  center  blocks,  directly  beneath  the 
vessel's  keel. 

5 .  Pier  -  A  column  built  of  blocks  that  extends  from  the 
ship  hull  to  dock  floor. 

6.  Side  Blocks  -  The  blocks  located  to  the  right  and  left  of 
the  keel  blocks. 

7.  Ton  -  A  long  ton,  2240  lbs. 

8.  g-  Acceleration  of  gravity,  32.2  ft/sec2. 


-17- 


SIDE   PIER 


KEEL    PIER 


FIGURE  2.1 
Vessel-Drydock  System 
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3.0   VESSEL-DRYDOCK  SYSTEM  FAILURE  MECHANISMS 

There  are  four  failure  modes  that  the  vessel-drydock 
system  can  exhibit:   block  crushing,  block  sliding,  block 
overturning  and  block  liftoff  (i.e.  the  liftoff  of  vessel  from 
a  row  of  blocks) . 

3 . 1  Block  Crushing6 

A  block  can  support  compressive  stresses  linearly  until 
the  blocks'  proportional  limit  is  exceeded.   The  proportional 
limit  is  a  material  property  of  the  block  and  is  the  maximum 
compressive  stress  at  which  stress  is  still  linearly 
proportional  to  strain.   If  the  compressive  proportional  limit 
is  exceeded,  the  block  is  considered  to  have  failed. 

3.2   Block  Sliding6 

Due  to  friction  at  the  interfaces  between  the  blocks  and 
the  vessel  and  between  blocks  with  other  blocks,  drydock 
blocks  will  have  the  tendency  to  resist  sliding  when  subjected 
to  vertical  and  horizontal  loads.   This  is  true  when  in  the 
absence  of  mechanical  fasteners. 

First  consider  the  blocks  that  form  the  keel  pier. 
Figure  3 . 1  shows  that  the  resistance  of  the  keel  pier  to 
sliding  is 

H  =  ylv    (block-block  interface) 
and  H  =  y 2V    (vessel-block  interface) 
where 

H  =  horizontal  resisting  force 
y i  =  coefficient  of  friction,  block-block  interface 
V-2   =   coefficient  of  friction,  block-vessel  interface 

V  =  vertical  load  on  the  blocks. 
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Failure  of  the  keel  pier  due  to  sliding  will  occur  when 
the  ratio  of  horizontal  force  to  the  vertical  force  is  greater 
than  the  lowest  coefficient  of  friction  found  in  the  block.  In 
other  words,  sliding  occurs  when 
|  >(v1   or  y2) • 

As  with  the  keel  piers,  the  resistance  of  the  side  piers 
to  sliding  will  depend  on  the  interface  which  offers  the  least 
resistance,  the  block-vessel  interface  or  the  block-block 
interface.   However,  this  case  is  more  complicated  due  to  the 
complex  geometry  of  the  side  pier  as  shown  in  Figure  3.2. 

The  block-vessel  interface  of  the  side  pier  is  examined 
and  shown  in  Figure  3.3.   An  arbitrary  force  F  at  some  angle 
from  vertical  is  applied  to  the  face  for  the  side  block  cap. 
At  this  interface,  the  horizontal  and  vertical  block  force 
reactions  are: 

Normal  Force  =  H  sin$  +  V  cos^ 
Tangential  Force  =  H  coscj)  -  V  sin<J> 
where  H  =  Horizontal  Reaction 

V  =  Vertical  Reaction 
$  =  Block  Inclination  Angle. 

Using  the  formula  for  horizontal  sliding  resistance,  H  =  y 
V,  previously  introduced  in  the  keel  pier  case  and  the 
above  equations  yields 

H  coscf)'-  V  sincp  =      y2    (H  sine})  +  V  cos<f>)  . 

Rearranging  the  above  equation  gives 

H         y?    cosi    +    sin6 


V        cosd>    -    u~    sin< 


or 


V  2 
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where  ^         p^    cos<$>   +   sin<J> 


y2  cos<t>    -    y2sin<P 

For  the  block-block  interface  of  the  side  pier,  sliding 
resistance  is  found  by  taking  the  vertical  and  horizontal 
components  of  the  applied  force,  F,  and  using  them  in  the 
formula  H  =  yj.  v*   Like  in  the  keel  pier,  failure  of  the  side 
pier  due  to  sliding  will  occur  when  the  following  is  satisfied 
|  >  (  yx  or  y2  ' )  . 

There  is  a  critical  angle  a  where  no  slippage  occurs. 
Referring  again  to  Figure  3.3,  this  angle  can  be  calculated  by 
comparing  the  maximum  tangential  force  F^  max  the  interface 
can  have  without  slipping  and  the  corresponding  normal  force 
Fn.   The  following  relationship  holds 

Ft  max  =  ^Fn 

Finally, 

a  =  arc  tan  y. 
Since  the  applied  force  F  acts  in  a  straight  line  through  all 
blocks  to  the  ground,  the  angle  oC  for  each  interface  must  be 
calculated.   If  force  F  is  applied  outside  of  these  angles, 
slippage  will  occur. 

3.3   Block  Overturnincf11 

The  third  failure  mechanism  is  the  overturning  of  a  block 
due  to  an  applied  force.   The  line  of  action  of  this  force 
must  fall  within  the  middle  one-third  of  the  base  of  the 
block,  as  shown  in  Figure  3.4,  or  the  block  will  tip  over. 

In  order  to  show  the  region  that  the  block  will  remain 
upright,  the  inclined  force  F  must  be  broken  into  a  transverse 
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component  Fn  and  a  vertical  component  Fv.   Superposing  the 

bending  and  axial  stresses  caused  by  F^  and  Fv  gives  the 

resultant  stress  at  any  point  in  the  block  as 

a    =    --J-   F   -  — - 
ax    lz   h    A 

where         L  =  block  height 

A  =  base  area,  here  bh 
Iz  =  base  moment  of  intertia,  here  ^^3/12. 
The  minimum  stress  point  on  the  bottom  of  the  block 
occurs  at  point  M  (Figure  3.4),  a  distance  -h/2  from  the  z 

axis.   Hence, 

r„  *     -Ik.  p Z 

lax'm   bh    h   bh   " 

When  axm  at  M  is  negative,  compressive  stress  is  present  and 
the  block  is  stable.   When  oxm  at  M  is  positive,  tensile 
stress  is  present  and  the  block  will  overturn  since  no 
fasteners  hold  block  to  the  ground  to  develop  tensile  stress. 
The  limit  of  block  stability  is  when 

(°x)m  =  0. 
Thus,  the  following  equation  holds 

F   =  A  F  . 
h    6L   v 

This  condition  exists  when  force  F  is  applied  at  an  angle  such 
that  the  line  of  force  lies  within  the  one-third  of  the  base 
of  the  block.   If  the  line  of  force  is  outside  this  region, 
the  block  will  tip. 

3.4   Vessel  Liftoff14 

The  fourth  and  final  failure  mode  occurs  when  the  vessel 
breaks  contact  with  either  the  side  or  keel  piers.   This 
failure  occurs  when  the  dynamic  deflection  of  a  row  of  blocks 
is  equal  to,  or  exceeds  the  average  static  deflection  of  the 
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blocks.   The  static  deflection  of  blocks  is  given  by 

«s      w 

s   2K    +  K, 

sv     kv 

where         W  =  submarine  weight 

Ksv  =  side  pier  vertical  stiffness 
Kkv  =  keel  pier  vertical  stiffness  « 

The  dynamic  deflection  can  be  due  to  vertical  displacement  y, 
rotational  displacement  v,    or  a  combination  of  both. 
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FIGURE  3.1 
Keel  Pier  Forces 


FIGURE  3.2 
Side  Pier  Forces,  General 
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FIGURE  3.3 
Side  Pier  Forces,  Block-Vessel  Interface 
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FIGURE  3.4 
Block  Tipping  Forces 
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4.0   MATHEMATICAL  MODELS  OF  VESSEL-DRYDOCK  SYSTEM 

This  portion  of  the  thesis  discusses  the  modelling  of  the 
drydocked  submarine  subjected  to  a  seismic  load.   Models  will 
include  a  quasi-static  method,  one,  two  and  three  degree  of 
freedom  mathematical  models.   The  models  developed  in  this 
section  will  be  analyzed  in  Section  6  and  7  to  predict  the 
vessel  seismic  response. 

4.1  The  Quasi-static  Method  of  Modelling  Seismic  Response 

4.1.1  Approach5 » 6 

Current  U.S.  Navy  design  method  for  submarine's  seismic 
response  will  be  examined  in  this  section.  The  submarine  is 
described  as  a  rigid  cylindrical  body  with  its  weight  evenly 
dispersed  longitudinally. 

The  quasi-static  force  method  replaces  the  earthquake 
motions  by  a  force  corresponding  to  the  vessel  mass  time 
0.2  g.   This  force  is  horizontally  applied  to  the  submarine's 
center  of  gravity  in  the  tranverse  direction.   The  drydock 
blocking  system  is  determined  by  the  evaluation  of  the  loads 
generated  by  the  quasi-static  force. 

4.1.2  Force  Equations 

The  application  of  the  static  force  to  the  center  of 
gravity  of  the  submarine  in  drydock  is  represented  in  Figure 
4.1.   The  seismic  overturning  moment,  Ms,  is  defined  as 

Ms  -  (A/9)  (a)  (KG)  (2240)  ft-lb 
where 

A  =  vessel  displacement  in  long  tons 

g  =  acceleration  of  gravity 

a  =  vessel's  center  of  gravity  acceleration. 
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Current  U.S.  Navy  design  practice  states  that  the  acceleration 

of  the  submarine's  center  of  gravity,  a,  should  be  set  equal 

to  0.2  g.   Hence 

Ms  =  448(A) (KG)  ft-lb 

After  determining  the  seismic  overturning  moment,  the 

number  of  sideblocks  required  on  one  side  of  the  vessel  is 

M 

N  =    S 


A  S  L 
P 


where        N  =  #  of  side  blocks  required  in  the  row 

A  =  contact  area  of  a  block  in  square  inches 

Sp  =  proportional  limit  of  the  block  cap  in  lb/in2 

L  =  distance  between  centerline  of  keel  and  side 
blocks,  in  feet,  as  shown  in  Figure  4.1. 

4.1.3   Vessel  Response  in  the  Quasi-static  Load 
Method 

The  application  of  a  force  applied  at  a  vessel's  center 
of  gravity  will  generate  reactionary  side  block  forces  in  the 
vertical  direction  to  oppose  the  seismic  moment,  as  shown  in 
Figure  4.2.   Summation  of  moments  about  the  keel  yields  the 
equation 

EM  =  Ms  -  L(Fsr  +  Fmr)  +  L(Fsl  -  Fml)  =0         (4.1.1) 


where 

Fs  =  static  pier  forces  (left  and  right) 

Fm  =  pier  forces  due  to  the  applied  moment  (left  and  right). 
Due  to  symmetry  in  the  vessel-drydock  system,     the 
applied  moment  resistive  forces  in  the  side  blocks  are  equal 
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in  magnitude  but  opposite  in  direction.   Thus,  Equation  4.1.1 
simplifies  to  the  form 

Ms  =  2L(Fm)  (4.1.2) 

where  Fm  =  Fml  =  Fmr. 

The  force  Fm  is  equal  to  the  side  pier  compression 
displacement,  6m,  times  its  vertical  stiffness,  Kvs.   This 
compression  displacement  can  be  expressed  as 

<5m  =  L  sin(6)  .  (4.1.3) 

where  6  =  angle  or  vessel  rotation  about  its  keel. 

Assuming  the  rotational  angle  is  small,  Equations  4.2.1  and 

4.1.3  are  combined  to  yield 

M 


2 
2L   Ksv  (4.1.4) 

4.2   One  Degree  of  Freedom;   Rotational  Response 
4.2.1  Approach 

A  one  degree  of  freedom  mathematical  model  of  a  submarine 
in  drydock  is  developed  in  this  section.   The  submarine  is 
identified  as  a  rigid  cylindrical  body  with  an  even 
longitudinal  weight  distribution.   The  system  damping  and 
stiffness  are  provided  by  the  drydock  blocking  arrangement. 
System  excitation  is  due  to  seismic  ground  accelerations. 
Other  assumptions  are 


The  keel  and  side  piers  remain  vertical  during  the 
ground  motions, 

no  slippage  between  the  block  cradle  and  drydock,  and 

the  keel  and  side  piers  bases  accelerate  at  the  same 
rate  as  the  drydock  (ground) . 

The  idealized  model  of  the  vessel-drydock  system  is  shown 
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in  Figure  4.3.   In  this  system,  vessel  motion  is  restricted  to 
rotation  about  the  keel.   The  side  piers  are  modelled  as 
vertical  due  to  typically  small  angles  or  inclination  of  the 
tops  of  the  side  piers.   Only  vertical  displacements  in  the 
side  piers  are  considered. 

The  one  degree  of  freedom  mathematical  model  predicts  and 
analyzes  vessel  rotational  response  about  its  keel  caused  by 
seismic  ground  accelerations  normal  to  the  submarine's 
longitudinal  axis.   The  resultant  angle  or  rotation,  expressed 
in  radians,  and  side  pier  forces,  generated  by  the  rotation, 
are  compared  to  applicable  pier  failure  criteria  found  in 
Section  3,  to  see  if  the  blocking  arrangment  will  remain 
intacted. 

4.2.2   One  Degree  of  Freedom  Equation  of  Motion  Model 
For  the  one  degree  of  freedom  case,  the  keel  is  the 
origin  with  clockwise  rotation,  6 ,  positive.   The  simplified 
vessel-block  system  showing  coordinate  6  at  rest  and  when 
excited  are  shown  in  Figures  4.3  and  4.4  respectively.   The 
moment  equilibrium  equation  about  the  keel  (origin)  is: 
ZMk  =  1^6  +  MKGXg  -  MKGyg  sin6 

=   B/2(Flsv   -    Fldv)    -    B/2(F2sv   +   F2dv)    +   WKG    sine   (4.2.1) 

where  Ik  =  Mass  moment  of  inertia  of  vessel  about  the  keel 

M  =  Vessel  Mass 
FlsvF2sv  =  static  side  block  forces  in  the  vertical  direction 

(Figure  4.4) 
FldvF2dv  =  dynamic  side  block  forces  in  the  vertical 
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direction  (Figure  4.4) 
W  =  vessel  weight  (Mg) 
and  XqrYq  =   seismic  ground  accelerations  in  their  respective 

directions 
This  equation  can  be  simplified  by  the  following: 
1.)   The  left  and  right  side  block  rows  have  the  same 

vertical  spring  constant,  Ksv/  have  the  same  static 
deflection  and  equal  but  opposite  dynamic  deflection  due 
to  symmetry,  i.e. 
Flsv  =  F2sv  =  Fsv  (4.2.2) 

and   |Fldv|=  |  F2dvl  =  Fdv  (4.2.3) 

2.)   The  dynamic  side  block  force,  Fdv,  as  seen  from  Figure 
4.5  equals  to  a  force  due  to  the  modelled  spring 
displacement  plus  a  dissipative  force  due  to  system 
damping.  In  other  words, 
Fdv (spring)  =  (B/2)KSV  sine 
and  Fdv (dissipative)  =  C(B/2) 6 
or  when  combined, 

Fdv  "  B/2  (Ksv  sin0  +  ce)  (4.2.4) 

where  Ksv  =  vertical  spring  constant  of  the  side  blocks 

C  =  vertical  damping  coefficient  of  the  side  blocks 
Substituting  Equations  4.2.2,  4.2.3,  and  4.2.4  into 
Equation  4.2.1  and  rearranging  components  yields 

•  •  • 

Ik6  +  (B2/2)C6   +  [(B2/2)KSV  -  W  KG]  sin  6 

=  -MKG  xg  +  MKG  yg  sin  6  (4.2.5) 

An  additional  simplification  is  that  since  Sis  very  small, 

sin  e  ~  e. 
The  final  simplification  is  that  damping  is  expressed  as  a 
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fraction  of  critical  damping.  To  do  this,  allow 

(B2/2)C  -  K       Ccr  =  2C  IfcOJn 
where  £  =  viscous  damping  factor(B2c/2Ccr) 

Ccr  =  critical  damping  coefficient  =  2Ikwn 

00 n  =  system  undamped  natural  frequency. 
Now,  Equation  4.2.5  simplified  to 

Ik6  +  2  K    IK%e  +  [(B2/2)KSV  -  WKG]6 

=  -MKG  xg  +  MKG  y g  6  (4.2.6) 

The  non-linear  term,  MKG  >;g9,  found  in  tne  equation  of 
motion  will  be  removed  at  this  time  to  linearize  Equation 
4.2.6  since  6  is  assumed  to  be  small.   The  effects  of  the 
non-linear  term  will  be  evaluated  in  Section  6  and  7. 

The  final  form  of  the  linearized  one  degree  of  freedom 
equation  of  motion  is 

•  •  • 

mjQ      +   C^e   +  k^   =  -m2  xg  (4.2.7) 

where  m^  =  Ik 

<=1  =  2?  Ik^n 
kx  =  (B2/2)KSV  -  W  KG 
m2  =  M  KG. 
In  matrix  form, 

o 


"o      m  " 

0 

• 

+ 

m1    Cl_ 

e 

-  - 

m,   o 


L"m2 


a  J 


To  simplify  the  explanation,  the  matrices  will  be  redefined  as 

[A]{yJ  +  [B]{y}  =  {E(t)}  (4.2.8) 

where 


[A]  = 


o   m. 


Lml  Cl 


[B]  = 


-ra. 


1-1 
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{y}  = 


(E(t)}  = 


-rru   x 
2   g 


Equations  4.2.6  and  4.2.8  will  be  solved  in  Sections  6  and  7. 

4.3   Two  Degree  of  Freedom;   Rotational  and  Horizontal 
Translation  Response 


4.3.1  Approach 

This  section  develops  a  two  degree  of  freedom 
mathematical  model  of  a  drydocked  submarine.   As  in  the  case 
of  the  one  degree  of  freedom  (Section  4.2.1),  the  model 
assumptions  still  hold,  i.e.  a  rigid  cylindrical  body,  keel 
and  side  piers  remain  vertical,  no  slippage  between  the  block 
cradle  and  drydock,  and  the  bases  of  the  keel  and  side  piers 
accelerate  at  the  same  rate  as  the  drydock  (ground) . 

The  idealized  model  of  the  vessel-drydock  system  is  shown 
in  Figure  4.6.   In  this  system,  vessel  motion  is  restricted 
to  rotation  about  the  keel,  and  in  horizontal  transverse 
translation  relative  to  the  keel  and  side  pier  supports.   The 
approach  for  determining  vessel-drydock  failue  is  identical  to 
the  one  degree  of  freedom  case  except  for  the  addition  of 
sliding  and  tipping  failure  modes. 

The  two  degree  of  freedom  mathematical  model  predicts  and 
analyzes  submarine  rotational  and  transitional  response  caused 
by  seismic  ground  accelerations  normal  to  the  vessel's 
longitudinal  axis.   The  rotational  response  is  in  radians,  the 
translational  response  in  inches,  and  both  responses  are 
expressed  as  a  function  of  time.   The  keel  and  side  pier 
forces  are  compared  to  applicable  pier  failure  criteria  found 
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in  Section  3 . 

4.3.2   Two  Degree  of  Freedom  Equation  of  Motion  Model 
As  in  the  one  degree  of  freedom  case,  the  keel  is  the 
origin  for  the  two  degree  of  freedom  system  with  clockwise 
rotation,  8  ,  is  positive.   The  linear  relative  translation 
coordinate,  x,  is  defined  with  respect  to  vessel  and  drydock 
translation  relative  to  the  ground.   Defining  u  as  the 
position  of  the  submarine  keel  relative  to  the  ground,  then 
the  following  relationships  hold: 

x  =  u  -  xg 

U   =   X   +   Xg 

and       u=x+Xg.  (4.3.1) 

The  simplified  vessel-block  system  showing  coordinates  (6  and 
x)  at  rest  and  when  excited  are  shown  in  Figures  4.6  and  4.7 
respectively. 

A  balance  of  forces  in  the  x  direction,  as  shown  in 
Figure  4.7,  yields  the  first  equation  of  motion 

ZFX  =  Mu  +  MKG6  =  (Flsh  -  Fldn)  -  (F2sh  +  F2dh) 

"  (F3sh  +  F3dh)    (4.3.2) 

where  Fisn,  F2sh,  F3sh  =  static  block  forces  in  the  horizontal 

direction  at  their  respective 
position 

Fldh'  F2dh'  F3dh  =  dynamic  block  forces  in  the 

horizontal  direction  at  their 
respective  position. 

The  above  equation  can  be  simplified  by  the  following: 
1)   The  left  and  right  side  block  rows  have  the  same 
horizontal  spring  constant,  Ksn,  have  equal  but  opposite 
static  horizontal  deflection  and  the  same  dynamic  horizontal 
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deflection  due  to  symmetry,  i.e., 

Flsh  +  F2sh  =  0 

Fldh  =  F2dh  •  (4.3.3) 

Also,  the  keel  block  rows  do  not  experience  any  net  static 
horizontal  deflection  since  the  inclination  angle  of  the  block 
cap  is  zero.   Hence 

F3sh  =  0.  (4.3.4) 

2)   The  dynamic  block  forces,  F^h  as  seen  from  Figure  4.7 
equals  to  a  force  due  to  the  modelled  springs  plus  a 
dissipative  force  due  to  system  damping.   In  other  words 

Fdh  (spring)  =  Fldh  (spring)  +  F2dh (spring)  +  F3dh  (spring) 
=  Kshx  +  Ksh  x  +  K^x 

Fdh  (spring)  =  (2Ksh  +  Kkh)x 
Likewise 

Fdh  (dissipative)  =  cxx 
so 

Fdh  =  (2Khs  +  Khk)x  +  cx*  (4.3.5) 

where  Ksn  =  side  pier  stiffness 
Kkh  =  keel  pier  stiffness 

Cx  =  system  horizontal  damping  coefficient. 
Substituting  Equations  4.3.1,  4.3.3,  4.3.4,  and  4.3.5 
into  Equation  4.3.2  and  rearranging  components  yields 

m^  x  +  m12  e  +  cl  x  +  kll  x  =  ~mll  xg  (4.3.6) 

where  m^  =  M 

m12  =  M  KG 

c^  =  Cx  =  system  horizontal  damping  coefficient 
kn  =  (2Ksn  +  K)^)  =  system  horizontal  stiffness 
Summing  the  moments  about  the  origin  in  a  similar  fashion 
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as  the  previous  degree  of  freedom  case,  as  shown  in  Figure 
4.7,  yields  the  second  equation  of  motion, 
IMK  =  1^6  +  MKG  u  -  MKG  yg  sine 

=  B/2(Flsv  -  Fldv)  -  B/2(F2sv  +  F2dv)  +  W*G  sin9  (4.3.7) 

The  second  equation  can  be  reduced  as  it  was  in  the  one  degree 
of  freedom  case  to 

Ike  +  MKG  x  +  Ce6    +  (B2/2  Ksv  -  WKG) 6 

=  MKG  Xg  +  MKG  yg6 .  (4.3.8) 

Once  again,  the  non-linear  term,  MKG  yg  6,  found  in 
Equation  4.3.8  will  be  removed  leaving  a  linearized  equation 
since  9  is  assumed  to  be  small.  The  effects  of  the  nonlinear 
term  will  be  evaluated  in  Sections  6  and  7. 

Now  the  system  of  equations  for  the  two  degree  of  freedom 
vessel-drydock  system  is  as  follows, 

mj!  x  +  m12e   +  en  x  +  kn  x  =  -i^  xg  (4.3.9a) 

and 

m226   +  m21  x   +    c126   +  k226  =  "m21  xg  (4.3.9b) 

where   m^  =  M 

m21  =  ^12  =  M^G 

m22  =  Ik 
c±   =  system  horizontal  damping  coefficient 

c2  =  system  rotational  damping  coefficient 

kll  =  2Ksh  +  Kkh 
k22  =  B2/2  Ksv  -  WKG. 
These  two  equations  of  motion  are  coupled  in  the  mass  times 
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acceleration  term,  which  is  known  as  inertial  coupling.   Also 
the  equations  are  coupled  in  the  damping  term,  because  the 
system's  two  natural  modes  have  translation  and  rotation 
involved  in  them.   This  model  uses  5  percent  of  the  critical 
damping  coefficient  as  system  damping.   To  evaluate  damping,  a 
modal  analysis  for  the  system  is  performed.  This  analysis  is 
performed  in  Appendix  3 .  Substituting  damping  coefficients 
into  Equations  4.3.9a  and  4.3.9b,  the  equations  of  motion  can 
be  expressed  in  matrix  notation  as  Equation  4.3.10: 

[M]  {y_i}  +  [C]  {y_!_}  +  [K]  {yj_}  =  (E'(t)}  (4.3.10) 
where 


{y1}  =  Response  Vector  = < A 
[M]   =  Mass  Matrix  = 


mll   m12 


Lm21 


m 


22J 


[C]   =  Damping  Matrix  = 


c,  ,   c,  ' 
li    lz 


Lc21   C22J 


[K]   =  Stiffness  Matrix  = 


kll   ° 


.o 


■2  2-1 


{E'(t)}  =  Seismic  Forcing  Vector  = 


-m,  ,  x 
li   g 

-m0,  x 
21   g 


To  simplify  the  explanation,  the  matrices  will  be  redefined  as 

[A]  {y}  +  [B]  {y}  =  {E(t)}  (4.3.11) 

where 


[A]  = 


o   [M] 
[M]   [C] 
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[B] 


{y} 


-[M] 
0 

{y'} 
{y'l 


o 

[K] 


{E(t)}  = 


O 


>{e-  (t)>;  " 


Equations  4.3.6,  4.3.8,  and  4.3.11  will  be  solved  in  Sections 

6  and  7. 

4.4   Three  Degree  of  Freedom:   Rotational,  Horizontal, 
and  Vertical  Translation  Response 

4.4.1  Approach 

This  section  develops  a  three  degree  of  freedom 
mathematical  model  of  a  drydocked  submarine.   The  idealized 
model  of  the  vessel-drydock  system  is  shown  in  Figure  4.8.   In 
this  system,  vessel  motion  is  restricted  to  rotation  about  the 
keel,  and  in  horizontal  transverse  and  vertical  translations 
relative  to  the  keel  and  side  pier  supports.   The  approach  for 
determining  vessel-drydock  failure  is  identical  to  the  two 
degree  of  freedom  case. 

The  three  degree  of  freedom  mathematical  model  predicts 
and  analyzes  submarine  rotational,  horizontal  and  vertical 
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transitional  responses  caused  by  seismic  ground  accelerations 
normal  to  the  vessel's  longitudinal  axis.   The  rotational 
response  is  in  radians,  the  translational  response  in  inches, 
and  all  three  responses  are  expressed  as  a  function  of  time. 
The  keel  and  side  pier  forces  generated  by  the  responses  are 
compared  to  applicable  pier  failure  criteria  found  in  section 
3  to  check  blocking  stability. 

4.4.2   Three  Degree  of  Freedom  Equation  of  Motion  Model 
As  in  the  two  degree  of  freedom  case,  the  keel  is  the 
origin  for  the  three  degree  of  freedom  system  with  absolute 
rotation,  6,  and  translation  coordinate,  x,  is  the  horizontal 
vessel  motion  relative  to  the  drydock.   The  linear  relative 
translation  coordinate,  y,  is  defined  with  respect  to  vessel 
and  drydock  translation  relative  to  the  ground.   This 
coordinate  y  describes  the  vertical  motion  of  the  vessel 
hull  relative  to  the  drydock  and  the  drydock  blocks.   Defining 
v  as  the  position  of  the  submarine  keel  relative  to  the 
ground,  then  the  following  relationships  hold; 

y  =  v  -  yg 

v  =  y  +  yg 
and  v=y+yg  (4.4.1) 

Likewise,  u  =  x  +  xg  (4.4.2) 

The  simplified  vessel-block  system  showing  coordinates  (x,  y 
and  6)  at  rest  and  when  excited  are  shown  in  Figure  4.8  and 
4 . 9  respectively . 

As  developed  in  Section  4.3.2,  the  first  and  third 
equations  of  motion  in  the  three  degree  of  freedom  model  are 

Mx  +   MKGO   +   Cxx  +    (2Ksh  +   K}ch)x  =   -Mxg  (4.4.3) 
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and 

Ik'e  +  MKG  x  -  MKGye  +  C© 

+  [(B2/2)KSV  -  WKG]6  =  -MKG  xg  +  MKG  yg6  .     (4.4.4) 

A  balance  of  forces  in  the  y  direction,  as  shown  in 
Figure  4.9,  yields  the  second  equation  of  motion 

lFy  =  Mv  =  (Flsv  -  Fldv)  +  (F2sv  -  F2dv) 

+  (F3sv  "  F3dv)  "  w-  (4.4.5) 

The  above  equations  can  be  simplified  by  the  following: 

1)  The  vessel  weight,  W,  must  be  equal  to  the  summation 
of  static  block  forces  so  that  static  equilibrium  holds  true, 
i.e. 

Flsv  +  F2sv  +  F3sv  =  w  (4.4.6) 

2)  The  dynamic  block  forces,  Fdv,  as  seen  from  Figure 
4.9  equals  to  a  force  due  to  the  modelled  springs  plus  a 
dissipative  force  due  to  system  damping.   These  dynamic  block 
forces  are  only  a  function  of  vertical  displacement,  y,  and 
velocity,  y.   The  formulation  of  these  forces  is  identical  to 
the  procedure  used  in  Section  4.3.2  and  will  not  be  repeated. 
The  final  form  of  the  dynamic  block  forces,  F^,  is 

Fdv  =  (2Ksv  +  Kkv)y  +  Cy  y  (4.4.7) 

where  Ksv  =  side  pier  vertical  stiffness 
Kkv  =  keel  pier  vertical  stiffness 
Cy  =  system  vertical  damping  coefficient. 
Substituting  Equations  4.4.1,  4.4.2,  4.4.6  and  4.4.7  into 
Equation  4.4.5  and  rearranging  components  yields 

My  +  Cyy  +  (2KSV  +  Kkv)y  =  -Myg.  (4.4.8) 

The  above  equation  along  with  Equations  4.4.3  and  4.4.4 
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formulate  the  three  equations  of  motion  for  the  three  degree 

of  freedom  model . 

The  non-linear  terms,  MKGyQ  and  MKGyge,  found  in 

Equation  4.4.4  will  be  removed  at  this  time  leaving  a 

linearized  equation  since  6  is  assumed  to  be  small.   The 

effects  of  the  non-linear  terms  will  be  evaluated  in  Section  6 

and  7. 

Now  the  system  of  equations  for  the  three  degree  of 

freedom  is  as  follows; 

rein   x  +  m13  9  +  cl  *  +  kll  x  =  ~mll  xg  (4.4.9a) 

m22  Y  +  c2  y  +  k22  y  =  _m22  yg  (4.4.9b) 

m33  e  +  m31  x  +  c3  6  +  k33 6  =  -m31  Xg         (4.4.9c) 

where  i^  =  m22  =  M 

m13  =  m31  =  MKG 

m33  "  xk 

ci  =  system  horizontal  damping  coefficient 

c2   =  system  vertical  damping  coefficient 

c3  =  system  rotational  damping  coefficient 

kll  =  2Ksh  +  Kkh 

k22  =  2Ksv  +  Kkv 

k33  =  B2/2  Ksv  -  WKG. 

The  first  and  third  equations  of  motions,  Equations  4.4.11a 

and  4.4.11c,  respectively,  are  coupled  in  the  mass  times 

acceleration  term,  which  is  known  as  inertial  coupling.   Also, 

these  equations  are  coupled  in  the  damping  term  because  two  of 

the  system's  three  natural  modes  have  horizontal  translation 

and  rotation  involved  in  them.   The  second  equation  of  motion, 

Equation  4.4.11b,  is  uncoupled  from  the  rest  of  equations* 
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Hence,  the  last  natural  mode  of  the  system  only  depends  on 
vertical  translation. 

In  order  to  introduce  damping,  this  model  uses  5  percent 
of  the  critical  damping  coefficient  as  system  damping.   To 
evaluate  damping,  a  modal  analysis  for  the  system  in  question 
must  be  performed.   This  analysis  is  performed  in  Appendix  3  . 
System  damping  coefficients  can  be  found  using  this  modal 
anlaysis  technique,  but  for  now,  the  coefficients  will  be  left 
in  terms  of  c^,  ci3/  c22'  c31  anc*  c33«   Substituting  these 
coefficients  into  Equations  4.4.11a,  4.4.11b,  and  4.4.11c,  the 
equations  of  motion  can  be  expressed  in  matrix  notation  as 
Equation  4.4.12: 

[M]  {y_i}  +  [C]  {yj.}  +  [K]  {yJJ  =  {E_!_(t)  }  (4.4.10) 
where 


{y'} 


Response  Vect 


Hi! 


[M] 


=  Mass  Matrix  = 


mll   ° 


m 


13 


m22   o 


31 


o 


m 


33 


[K] 


=  Stiffness  Matrix  = 


k,  ,   o     o 


o     k22   o 


{E'(t)}  =  Seismic  Forcing  Vector  =  j  — m-,  -. 


33 


x 


-m00  y 

22  J  g 

-m-.,  x 
31    g 


[C] 


=  Damping  Matrix  = 


Cll   °     C13 


c22   o 


c31   ° 


33 
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To  simplify  the  explanation  the  matrices  will  be  redefined  as 

[A]  {y}  +  [B]  {y}  =  (E(t)}  (4.4.11) 

where 


[A]  = 


[B] 


0 

[M] 

[M] 

[C] 

- 

-[M] 

0 

0 

[K] 

{y}  = 


{E(t)}  = 


11} 


\r\ 


{E' (t)} 


Equations  4.4.3,  4.4.4,  4.4.8,  and  4.4.11  will  be  solved  in 
Section  6  and  7. 
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FIGURE  4.1 
Vessel-Drydock  System  with  Applied  Quasi-  static  Force 


Fs  -F„ 


F"S+FM 


FIGURE  4.2 
Vessel-Drydock  System  Force  Diagram  with  Applied  Quasi-static 

Force 
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FIGURE  4.3 
Idealized  One  Degree  of  Freedom  Vessel-  Drydock  System  at  Rest 
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FIGURE  4.4 
Idealized  One  Degree  of  Freedom  Vessel-  Drydock  System, 

Excited 
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FIGURE  4.5 
Dynamic  Side  Block  Force 
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FIGURE  4.6 
Idealized  Two  Degree  of  Freedom  Vessel-Drydock  System  at  Rest 
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FIGURE  4.7 
Idealized  Two  Degree  of  Freedom  Vessel-Drydock  System, 

Excited 
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FIGURE    4.8 
Idealized  Three   Degree   of   Freedom  Vessel-Drydock  System  at 

Rest 
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FIGURE  4.9 
Idealized  Three  Degree  of  Freedom  Vessel-Drydock  System, 

Excited 
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5.0   SYSTEM  PARAMETERS 

In  order  to  study  the  response  of  a  vessel-drydock  system 
to  seismic  loading,  eleven  typical  configurations  have  been 
selected  for  implementation  into  the  one,  two  and  three 
degrees  of  freedom  models  developed  in  the  previous  section. 
This  section  explores  these  eleven  systems  along  with  the 
modeling  of  the  pier  blocks  for  stiffness. 

5.1  Vessel-Drydock  System  Parameters 

Eleven  drydock  system  configurations  have  been  chosen  for 

analysis  and  are  defined  in  Table  5.1.   The  vessel-drydock 

parameters,  corresponding  to  each  system,  were  taken  from  the 

appropriate  NAVSEA  drawings. 

TABLE  5.1 

Vessel-Drydock  Configurations 

Longitudinal     NAVSEA 
Block  Spring Drawing 

845-2006640 

845-2006640 

845-2006640 

845-2006640 

845-2006640 

845-4862749 
845-4862749 
845-4862749 
845-4403511 
845-2140554 
845-2140554 


Block 

stem 

Hull 

Type 

1 

616 

Composite 

2 

616 

Composite 

3 

616 

Timber 

4 

616 

Timber 

5 

616 

Timber  Side/ 
Concrete  Keel 

6 

726 

Composite 

7 

726 

Composite 

8 

726 

Composite 

9 

688 

Composite/Timber 

10 

637 

Composite/Timber 

11 

637 

Composite/Timber 

8 

ft 

16 

ft 

8 

ft 

16 

ft 

16 

ft 

8 

ft 

12 

ft 

16 

ft 

12 

ft 

12 

ft 

12 

ft 
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Table  5.2  lists  the  vessel  radii,  weight,  distance  from  center 
of  gravity  to  keel,  KG,  moments  of  inertia  about  vessel  center 
of  gravity,  Ic.g.  and  keel,  1^,  and  side  block  transverse 
spacing,  B. 

TABLE  5.2 

Key  Vessel  Parameters 

Hull 
SYSTEM  Radius   Weight    KG        Icg        Ik  B 

(inches)  (Kips)  (inches)  (K-in-sec2)  (K-in-sec2)  (inches) 


1-5      198     16,396    193       831,257      2,411,000     144 

6-8      252     37,656    223     3,097,535      7,949,000     180 

9       198     13,624    193       691,852      2,007,000     138 

10-11     190      9,529    174       445,592      1,193,000     138 

5.2   Block  Parameters11 

The  models  developed  in  Section  4  require  the  evaluation 
of  the  block's  vertical  and  horizontal  spring  stiffness 
constants,  Kv  and  Kn.   The  spring  stiffness  is  the  amount  of 
force  required  to  produce  one  unit  of  displacement  in  the 
respective  direction.   In  order  to  facilitate  the  use  of 
Hooke's  law  in  evaluating  spring  constants,  the  assumption  of 
homogeneous  isotropic  behavior  will  be  used  for  the  block 
materials:   wood  and  concrete.   Both  the  side  and  keel  piers 
must  be  analyzed  for  Kv  and  K^. 

5.2.1  Vertical  Stiffness  KV1:L 
Using  Hooke's  Law,  the  linear  relationship  between 
vertical  applied  force  and  deformation  can  be  written  as 

F  _  Eh 

A  "  H 
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where  F  =  vertical  force 

A  =  area  under  applied  force 
E  =  modulus  of  elasticity 
h  =  change  in  height 
H  =  original  height. 

Hence,  the  spring  constant  for  a  given  material  is 

F    EA 
v    h  '   H   » 

(5.1) 

For  the  eleven  drydock  configurations,  two  types  of 

standard  blocks  (composite  and  timber)  are  used.   A  standard 

composite  block  which  is  composed  of  a  softwood  cap,  a 

hardwood  middle  portion,  and  a  concrete  bottom  can  be  modelled 

as  a  series  of  three  springs  (Figure  5.1).   The  resultant 

block  spring  stiffness  is 

kv  block  =  [(lAcap)  +  (V^oak)  +  UAcon)]-1     (5-2) 
where 

kcap  =  softwood  cap  spring  constant 

koak  =  hardwood  middle  spring  constant 

kcon  =  concrete  bottom  spring  constant. 

All  three  constants  listed  above  are  calculated  from 
Equation  5.1.   A  standard  timber  block  is  composed  of  a  softwood 
cap  and  a  hardwood  body  (Figure  5.2).   Its  block  spring 
constant  is 

*v  block  =  [(lACap)  +  UAoak)]-1.  (5.3) 

The  total  vertical  stiffness  for  a  drydock  configuration 
can  be  computed  by  multiplying  kv  block  times  the  number  of 
blocks.   The  keel  pier  has  a  total  vetical  stiffness  of 
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Kkv  "  kv  block  *  nk 
and  for  the  side  pier 

Ksv  =  kv  block  *  ns 
whre  n^  and  ns  are  the  number  of  blocks  in  one  row  for  their 

respective  pier. 

The  values  of  total  vertical  stiffness  for  both  keel  and 
side  blocks  of  the  eleven  drydock  configurations  are  found  in 
Table  5.4.   The  following  parameters  are  used  in  the 
formulation  of  Table  5.4: 

A  =  cross-sectional  area  of  softwood  cap 

Ecap  =22.5  ksi 

Eoak  =  31.675  ksi 

Econ  =  2,000  ksi 

Hcap  =  4  inches 

Hoak  =  29  inches  for  composite  block 
=  56  inches  for  timber  block 

Hcon  =27  inches. 

5.2.2  Horizontal  Stiffness  K^11 

In  order  to  determine  the  horizontal  stiffness,  two  types 
of  deformation  must  be  looked  at.   They  are  the  block-cap 
displacements  due  to  bending  and  shear  deformations.   Modelling 
as  a  continuous  cantilever  beam  subjected  to  a  concentrated 
lateral  force  at  the  cap  surface  (Figure  5.3),  the  bending 
displacement  due  to  the  applied  force  P  of  a  composite  block  is 

P(H;,+H2+H3)    P(E1-E2)  (H2+H3) 
d_  =  =■=  ^ + 


B  3E,E2I 
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P(E2I-E3I3)H3 
+        3E2E3II3  (5>4) 

where  H^  =   height  of  concrete  (27  inches) 

H2  =  height  of  hardwood  (29  inches) 

H3  =  height  of  softwood  cap  (4  inches) 

Ej  =  modulus  of  elasticity  of  concrete 

E2  =  modulus  of  elasticity  of  hardwood 

E3  =  modulus  of  elasticity  of  softwood  cap 
I  =  moment  of  inertia  of  block's  cross-section 

I3  =  moment  of  inertia  of  cap's  cross-section. 
For  the  bending  displacement  of  a  timber  block,  Equation  5.4 
still  holds  true  with  the  alterations: 

K±   =  height  of  hardwood  (27  inches) 

E]_  =  modulus  of  elasticity  of  hardwood. 
In  shear,  deformation  can  be  determined  by  modelling  the 
composite  block  as  an  element  subjected  to  shear  stress  at  the 

top  (Figure  5.4).   The  shear  displacement  becomes 

(1+v         )PH,  2(l+v  JPHn 

-,      _   con         1  wood         2 

s  '       AE  AE„ 


2(1+Vwood)PH3 
AE3 


(5.5) 

where  A  =  area  of  cap's  cross  section 

vwood  =  Poisson's  ratio  for  wood  (0.30  is  used) 
vcon  =  Poisson's  ratio  for  concrete  (0.15  is  used). 

By  redefining  H^  and  E±,    Equation  5.5  yields  the  shear 

displacement  of  a  timber  block. 
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The  horizontal  spring  constant  of  any  one  block  can  be 
determined  through  the  Hooke's  Law  relation 

kh  =  P/(dB  +  ds)  (5.6) 

The  total  horizontal  stiffnesses  for  the  keel  and  side  piers 
are 

Kkh  "  kh  *  nk 
and  Ksh  =  kh  *  ns  . 

The  values  for  the  total  area  of  the  rows  of  blocks  for 

each  of  the  eleven  systems  are  listed  in  Table  5.3.   Also  listed 

are  the  corresponding  total  area  of  the  cap  blocks  for  a  row  of 

side  or  keel  blocks.   The  values  of  total  horizontal  and 

vertical  stiffnesses  of  the  keel  and  side  piers  for  the  eleven 

drydock  configurations  are  listed  in  Table  5.4. 
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TABLE  5.4 
Total  Keel  and  Side  Pier  Stiffness  in  Pounds  Per  Inch 


System 

Kkv 

KSV 

*kh 

Ksh 

1 

50089159 

7545899. 

14018741 

2612721 

2 

50089159 

3903051. 

14018741 

1351407. 

3 

28493109 

4292468. 

5161684. 

1282216. 

4 

28493109 

2220242. 

5161684. 

663215.3 

5 

50089159 

2220242. 

14018741. 

663215.3 

6 

98356895 

52105736 

24330635 

16104611 

7 

98356895 

34542004 

24330635 

10676090 

8 

98356895 

26345596 

24330635 

8142780. 

9 

21758374 

8239775. 

3941650. 

2683291. 

10 

16577809 

6765289. 

3072572. 

2139500. 

11 

16577809 

5464271. 

3072572. 

1728057. 
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SQFTVQQD  CAP 


HARDVQDD 


CAP 


FIGURE  5.1 
Standard  Composite  Block  Stiffness  Model 


SDFTVDDD  CAP 


CAP 


FIGURE  5.2 
Standard  Timber  Block  Stiffness  Model 
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FIGURE  5.3 
Standard  Block  in  Bending 


FIGURE  5.4 
Standard  Block  in  Shear 
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6.0   VESSEL  RESPONSE  TO  SEISMIC  EXCITATION 

This  portion  of  the  thesis  discusses  methods  of 
determining  the  seismic  response  of  the  vessel-drydock  system 
mathematical  models  developed  in  Section  4 .   The  numerical 
analysis  schemes  to  be  implemented  are  a  modal  analysis  method 
(Section  6.1)  which  looks  at  the  linearized  equations  of 
motion,  and  a  fourth  order  Runge-Kutta  numerical  method 
(Section  6.2)  which  analyzes  the  non-linear  aspect  and  also 
verifies  the  linear  solution.   Section  6.3  incorporates  a 
response  spectrum  analysis  to  solve  the  mathematical  models. 
Finally,  Section  6.4  generates  system  response  for  the  eleven 
vessel-drydock  configurations  using  the  quasi-static  force 
method. 

6. 1  Vessel  Response  Using  Modal  Analysis 
In  order  to  use  the  modal  analysis  method  described 
herein,  only  the  linearized  equations  of  motion  mathematical 
models  can  be  evaluated.   These  are  Equations  4.2.8,  4.3.11 
and  4.4.11  which  correspond  to  the  one,  two,  and  three  degree 
of  freedom  models  respectively.   Since  all  three  linearized 
equations  are  in  the  form 

CA]{y}  +  [B]{y}  =  (E(t)}, 
the  three  degree  of  freedom  model  will  be  analyzed.   The  other 
two  cases  are  just  a  reduced  version  of  the  third. 

To  solve  the  linearized  three  degree  of  freedom 
equations  of  motion  (4.4.11)  in  matrix  notation  and  obtain 
system  responses  as  a  function  of  time,  the  problem  is 
solved  as  an  eigenvalue  problem  with  a  procedure  described  in 
reference  [13].   As  a  first  step,  the  equations  are  decoupled. 
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To  accomplish  this,  first  consider  the  free  vibration  problem 
[A]{y}  +  [B]{y)  =  0  (6.1.1) 

Assuming  Equation  6.1.1  has  a  solution  in  the  form 

{y}  =  {^}  eAt  , 
Equation  6.1.1  can  be  expressed  as 

MA]{i|>)  +  [B]{i|j}  =  0.  (6.1.2) 

where   X  =  complex  eigenvalues 

{  i>}   =  modal  vectors. 
Solving  the  matrix  determinant  problem 

I  X[A]  +  [B] |  =0  (6.1.3) 

yields  the  eigenvalues.   In  the  three  degree  of  freedom  case, 
[A]  and  [B]  are  6x6  matrices  hence  the  solution  of  Equation 
6.1.3  will  yield  6  eigenvalues.   These  eigenvalues  will  be 
comprised  of  three  sets  of  complex  conjugates.   Now  the  modal 
vectors  can  be  found  by  substituting  the  eigenvalues  back  into 
Equation  6.1.2  and  solving.   The  modal  vectors  are  then 
combined  to  form  the  modal  matrix  [  ty]  where 

[*]  =  [(Wi  (iW2  (^3  U>4  {^}5  <iJ>)6]' 
The  equations  of  motion  can  then  be  decoupled  by  using  the 
modal  matrix,  i.e. 

U]T  [AH1>]  =  C^A^.] 
and 

C^]T  [B][i|>]  =  C— B_] 

where  [ — A-_]  and  [ — -B .]  are  diagonal  matrices. 

Once  the  modal  matrix  is  found,  the  forced  vibration 
problem  is  analyzed.   Assume  the  following 

(Y)  =  W    <z) 
and  {z}  =  m"1  {y}, 
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where  z  is  a  complex  coordinate  system.   Equation  4.4.1  is  now 
expressed  as 

[A]     [HO    {z}    +    [B]     [i|i]    {z}    =    (E(t)}. 
Using    [^]T,    the  above  equation  can  now  be   in  the   form 

r-A^]  {z}  +  [-—  B— ]  {z}  =  mT{E}  ={N(t)}. 
This  will  yield  six  uncoupled  equation  of  motion  in  the  form 

aii  Zi(t)  +  bii  Zi(t)  =  N^t),   i  =  1,6 
or   z'i(t)  -  vi  Zi(t)  =  (l/aii)  Ni(t)  (6.1.4) 

where  x\  =   -b^i/a-j^. 
The  solution  of  Equation  6.1.4  is 


Zi(t)  =  /  e  vi(t-T)  Ni(T)  dx,  i  =  1,6.      (6.1.5) 

11  o 
Once  Equation  6.1.5  is  solved,  the  {y}  matrix  is  found  by  the 

relationship 

{y.}  =  m  {z}. 
Recalling  from  Section  4.4.2  the  terms  of  the  three 
degree  of  freedom  equations  of  motion  {y}  matrix  is 

x 


{y.}  = 


Hence,  the  desired  vessel-drydock  system  response  can  be  found 
in  the  fourth,  fifth  and  sixth  terms  of  the  {y}  matrix,  i.e. 
x,  y  and  6. 
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6.2   Vessel  Response;   Fourth  Order  Runae-Kutta  Numerical 
Method 

A  fourth  order  Runge-Kutta  numerical  scheme  is  chosen  to 
analyze  the  non-linear  equations  of  motion  due  to  its 
simplicity  and  easy  use  of  computer  programming.   It  also  will 
be  used  to  verify  solutions  found  by  other  methods.   The 
non-1 inear  Equations  4.2.6;  4.3.6  and  4.3.8;  4.4.3.  4.4.4  and 
4.4.8  correspond  to  the  one,  two  and  three  degree  of 
freedom  models  respectively.  Since  all  three  systems  of 
equations  are  similar  in  form,  only  the  three  degree  of 
freedom  non-linear  model  will  be  developed. 

The  three  degree  of  freedom  non-linear  equations  of 
motion  are 

mu   x  +  HI136      +   Cn   x  +   c13   6    +  ki:L  x  =   -m^   xg 

m22  y  +  c22  y  +  k22  y  =  -m22  yg 

and 

•  •  • 

m33e   +  m31  x  -  m31  y6+  c33e   +  c31  x  +  k31e 

=  -m31  xg  +  m31  yg  6. 
The  coefficient  of  the  above  equations  can  be  found  in  Section 
4.4.2.   These  equations  can  be  rearranged  into 

l  +  l!l3S-Jlii.!l3:-^ix-x_  (6.1.2a) 


c  y 

c22  •    22 

y  = v y  -  y 

m22  *   m22       g 


(6.1.2b) 


and   6  +  !!l31  "    _!33  J  _  ^31  '   ,^33    m31C22  •  +  m31k22   x 
m33      m33     m33     V33    m33m22     m22m33 


»3i5 

m33   g 


(6.2.1c) 
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Equations  6.2.1a,  6.2.1b  and  6.2.1c  can  be  made  into  an 
equivalent  first  order  system  by  the  following  substitutions 

x  =  R  =  ei(t,x,R,T) 

y  =  S  =  f!(t,y,S) 

e  =  t  =  gi(t,e,y,R,s,T) 

m13  • 
R  +  ^   T  =  e2 (t,x,R,T) 

S  =  f2  (t,y,S) 

m31 
T  +  -^  R  =  g  (t,  6,y,R,S,T)  . 

m13 


The  Runge-Kutta  formulas  [8]  for  this  system  are 

fell   =  hRn  £11   =  h{e2(t,xn/Rn,Tn)} 

fell  =  hSn  l12   =  h<f2(t,yn/Sn)} 

fell  =  hTn  £13   =  h{g2(t,en,    yn/Rn.sn^Tn)} 

fe21   =  h(Rn  +2fn)  l2l   =  h{e2(t,xn  +   2fell^n  +  kll^n   +  Wll3)  ) 

1  1 

fe22    =   h(Sn   +2^12)  £22    =   h{f2(t/yn   +    2fe12/sn   +   Ull)  > 

1  1 

fe23    =  h(Tn  +2^13)  ^23    =  h{g2(t,6n   +   2fe13/yn  +j=>fe12' 

Rn  +   ^lli    Sn  +   lfc12,    Tn  +  lt13)  } 

fe31   =  h(Rn   +    2^21)^31   =  h{e2(t,xn  +   3*21**11  +  l*21f*ll  -^23)} 

1  1 

fe32   =  h(Sn  +    2^22)£32   =  h{f2(t,yn  +   2fe22>sn  +  i*22> > 

1  1 

fe33    =  h(Tn  +    2^23)^33    =  h{g2(t/6n  +   2*23^11  +  ^-fe22/Rn  +  ^21> 

1  1 

sn  +    2£22'Tn  +   2£23> > 
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fe41  =  h(Rn  +  l31) 
fe42  =  h(Sn  +  l32) 
fe43    =   MTn   +    ^3) 


and  lastly, 

xn+l  = 

xn  + 

1, 
6^ 

Yn+1  - 

yn  + 

1, 

Qn+l  = 

en  + 

1. 

Rn+1  = 

Rn  + 

i. 

sn+l  = 

sn  + 

i( 

Tn+1  = 

Tn  + 

1( 

lA1   =  h{e2(t,xn  +  fe31,  Rn  +  l31,Tn   +  ^33)} 

^42  =  h{f2(t,yn  +  k22l  Sn  +  ^32)} 

£43  =  h{g2(t/en  +  fe33>  Yn  +  fe32'Rn  +  £31» 

Sn  +  ^32/  Tn  +  £33)> 


fell  +  2fe21  +  2fe31  +  fe4l) 


+    2fc22  +  2'-32  +  R42) 
-13  +  2fc23  +  2fe33  +  fe43) 


lll   +    2^21  +  2^31  +  lA1) 


+    2l22    +    2l32    +    l42) 
+  2^2  3  +  2^33  +  £43) 


c 


12 


(' 


13 


where  h  is  the  time  step.   Hence,  the  desired  vessel-drydock 
system  response  can  be  found  and  system  failure  analyzed  by 
the  iteration  of  the  Runge-Kutta  scheme. 

6.3  Vessel  Response:  Response  Spectrum  Analysis 
The  degrees  of  freedom  mathematical  models  developed  in 
previous  sections  produce  a  time  history  of  the  seismic 
response  of  a  drydocked  vessel.   These  models  use  numerous 
computer  iterations  to  find  their  time  history  in  order  to 
evaluate  the  system's  maximum  seismic  response.   Another  way 
to  determine  the  maximum  seismic  response  of  a  drydocked 
vessel  is  the  response  spectrum  analysis  method.  Simplicity 
and  fast  results  are  characteristic  of  this  spectrum  method. 
However,  only  maximum  response  is  produced  and  no  time  history 
is  generated. 

The  response  spectra  are  graphs  of  the  maximum  seismic 
response  of  single  degree  of  freedom  systems  over  a  range  of 
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system  natural  frequencies.   These  graphs  are  generated  for 

particular  earthquake  and  can  contain  acceleration,  velocity, 

displacement,  or  combination  of  all  three  responses  for  the 

given  system.   Simplified  response  spectra  are  plotted  for 

earthquake  design  analysis.  The  El  Centro  earthquake  response 

spectra  graph,  shown  in  Figure  6.1,  will  be  used  to  confirm 

the  results  of  the  mathematical  models  in  this  thesis.   The 

simplified  spectrum  portion  of  this  graph  is  described  in 

reference  [2]  and  updated  in  reference  [7], 

The  simplified  El  Centro  earthquake  response  spectra  used 

described  the  maximum  relative  displacement  of  the  system 

shown  in  the  insert  of  Figure  6.1.   The  linear  equation  of 

motion  for  this  simple  system  is 

mu  +  cu  +  ku  =  -m  ysa  fa(t) 

where         u(t)  =  relative  displacement  of  the  mass  with 

respect  to  the  support 

and  ysa  =  maximum  support  acceleration. 

The  natural  frequency  of  the  system  is  given  by 

oj2  =  k/m 

f  =  W/27T   . 

Using  the  response  spectra  graph  and  the  particular  frequency, 
f,  of  the  system  of  interest,  the  maximum  value  of  the 
relative  displacement,  um,  is  determined.   In  equation  form, 
the  relative  displacements  of  a  single  degree  of  freedom 
system  with  5%  critical  damping  subjected  to  the  El  Centro 
earthquake  ground  motions  are: 
for  oj  <  2.24  rad/sec, 

umax  =  !-4  (Ysa)max  =  H-62  inches,       (6.3.1) 
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for   2.24  <  w  <  12.74, 

umax  =1.9  (Ysa)max/w  =  26.03/eo  inches,    (6.3.2) 

for  u  >  12.74, 

umax  =  2-6  (ySa)max/w2  =  331.53/0)2  inches.  (6.3.3) 
These  formulas  generated  the  simplified  response  spectrum 
shown  in  Figure  6.2. 

In  order  to  use  Figure  6.2  to  verify  the  one,  two  and 
three  degree  of  freedom  vessel-drydock  models,  some 
manipulation  must  be  performed.   This  is  because  the  linear 
equations  of  motion  of  the  models  do  not  exactly  match  those 
of  the  system  pictured  in  Figure  6.1.   Examining  the  one 
degree  of  freedom  linear  equation  of  motion  (Equation  4.2.7); 

Ik6  +  2Clkwn6  +[(B2/2)KSV  -  WKG]6  =  -M  KG  Xg 
or 

6  +  2?oj  n6  +  ([(B2/2)KSV  -  W  KG]/Ik)e  =  -  (M  KG/Ik)  xg. 
The  natural  frequency  for  the  above  equation  is  defined  as 

wn2  =  ((B2/2)KSV  -  W  KG)/Ik. 
Using  Figure  6.2,  this  frequency  gives  a  corresponding  umax. 
Hence,  the  maximum  rotational  response  of  the  system,   ,  is 
then  determined  by  the  relation 

emax  =  M  KG/!k  umax- 
Table  6.1  lists  the  maximum  rotational  response  for  the  eleven 

configurations  found  in  Section  5. 
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TABLE  6.1 

One  Degree  of  Freedom  Vessel-Drydock  System 
Maximum  Response  Using  Response  Spectrum  Analysis 


System  wn 


(rad/sec) 


umax 

frads. ) 

0. 

01594 

0. 

,02619 

0. 

02149 

0. 

03100 

0. 

03100 

0. 

,00698 

0. 

,00860 

0. 

,00987 

0. 

,01445 

0. 

,01299 

0. 

,01450 

1  5.580 

2  3.933 

3  4.141 

4  2.870 

5  2.870 

6  10.254 

7  8.327 

8  7.255 

9  6.147 

10  7.253 

11  6.498 


The  response  spectrum  analysis  to  the  two  and  three 
degree  of  freedom  systems  is  applied  by  means  of  the  modal 
method  described  in  reference  [2].   Concisely,  the  modal 
method  requires  that  the  system  natural  modes  be  determined. 
Then,  the  response  of  the  system  to  a  known  forcing  input  can 
be  developed  by  mode  superposition.   Treating  each  natural 
mode  as  a  single  degree  of  freedom  system  with  its  own  natural 
frequency,  this  method  is  good  for  simplifying  seismic 
response  analysis.   Since  each  natural  mode  acts  as  a  separate 
system,  the  maximum  response  may  be  determined.   Then,  a 
conservative  estimate  for  the  maximum  response  of  the  original 
system  can  be  made  by  adding  up  the  maximum  response  of  each 
mode. 

The  method  of  adding  maximum  modal  responses  gives  an 
upper  bound  for  the  maximum  system  response.   Since  there  is  a 
known  response  spectrum,  the  amount  in  which  each  natural  mode 
contributes  to  the  maximum  response  for  a  given  input  can  be 
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determined.   This  is  known  as  the  participation  factor  of  each 
given  mode.   A  full  description  and  derivation  of  partici- 
pation factors  including  implementation  techniques  can  be 
found  in  reference  [2],  and  will  not  be  presented  here.   The 
participation  factor  for  each  natural  mode  in  the  two  and 
three  degree  of  freedom  case  is  calculated  to  confirm  computer 
generated  results  of  maximum  responses  in  Appendix  3. 

6. 4   Vessel  Response:   Quasi-Static  Force  Analysis 
As  described  in  Section  4.1,  the  quasi-static  force 
method  replaces  the  earthquake  motions  by  a  force 
corresponding  to  the  vessel  mass  times  0.2g.   This  analysis 
only  allows  rotation  about  the  vessel's  keel  hence,  the  quasi- 
static  force  analysis  is  itself  a  one  degree  of  freedom  model 
with  system  response,  in  radians,  determined  by  Equation 
4.1.4,  i.e. 

B  =  Ms/L2  Ksv 


where      Ms  =  448  *  A  *  KG 

A  =  displacement  in  tons 
L  =  B/2. 
The  results  of  this  method  are  listed  in  Table  6.2. 
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TABLE  6.2 
System  Response  using  Quasi-Static  Force  Analysis 

SYSTEM  e,  IN  RADIANS 

1  0.00809 

2  0.01564 

3  0.01422 

4  0.02749 

5  0.02749 

6  0.00199 

7  0.00300 

8  0.04419 

9  0.00670 

10  0.00051 

11  0.00637 
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FIGURE  6.1 
Actual  El  Centro  Earthquake  Response  Spectrum 
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FIGURE  6.2 
Idealized  El  Centro  Earthquake  Response  Spectrum,  5%  Damping 
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7.0   COMPUTER  SOLUTIONS 

In  this  section  of  the  thesis,  the  computer  programs 
needed  to  implement  the  mathematical  models  and  the  associated 
seismic  response  are  developed.   Due  to  the  numerous 
iterations  that  is  required  for  the  modal  analysis  and 
Runge-Kutta  methods  to  effectively  determine  vessel-drydock 
response,  all  computer  work  is  performed  at  the  MIT  Joint 
Computer  Facility  (JCF) .   The  required  computer  programs  are 
written  in  Fortran  77. 

In  order  to  check  the  computer  generated  results,  the 
vessel-drydock  mathematical  models  are  subjected  to  a 
sinusoidal  earthquake  input.   The  advantage  being  that  closed 
form  solutions  of  responses  to  sinusoidal  inputs  can  be 
calculated.   Thus,  closed  form  response  can  be  compared  to  the 
computer  generated  responses. 

Upon  obtaining  the  correct  results  for  a  sinusoidal 
earthquake,  the  El  Centro  earthquake's  acceleration  time 
history  is  applied  to  the  vessel-drydock  system.   The  maximum 
value  of  these  computer  generated  results  is  compared  to  the 
maximum  response  for  the  system  predicted  by  the  response 
spectrum  method  found  in  Section  6.3. 

7 . 1  Modal  Analysis  Solution 

7.1.1  Computer  Program  Development 

As  described  in  Section  6.3,  the  program  approaches  the 
problem  by  first  assembling  matrices  A  and  B.   The  program 
then  uses  the  International  Mathematical  and  Statistical 
Library  (IMSL)  subroutine  EIGZF  resident  on  JCF  to  perform  the 

-73- 


system  free  vibration  eigenvalue  problem  on  matrices  A  and  B. 
The  output  of  subroutine  EIGZF  is  the  modal  matrix  [ip]  of  the 
vessel-drydock  system.   Reference  [4]  contains  further 
information  about  subroutine  EIGZF. 

The  equations  of  motion  can  now  be  de-coupled  by  using 
the  modal  matrix.   De-coupling  the  equations  of  motion  leaves 
equations  of  the  form 

aii  zi(t)  +  bii  Zi(t)  =  Nj.(t)  ,   i  =  1,  m 
where  m  is  equal  to  twice  the  number  of  the  particular  degree 
of  freedom  model,  i.e.  for  the  one  degree  of  freedom  ~  m  =  2. 
The  solution  to  the  above  equation  is 

t 

"bii/aii(t-T) 

zi(t)  =  er~     I      e  Ni(T)  dT'  i  =  l#»i 


-*/ 


o 
where  a^j_  and  bj_j_  represent  the  diagonal  elements  of  the 
diagonalized  A  and  B  matrices,  respectively.   The  integral 

t 


/ 


"bii/aii  (t"T> 
e 


N. (t)  dT 

i 


can  be  evaluated  using  Simpson's  rule,  which  is 

ZAi(t)  =  [(Aj_(t  -  2At)  +  Ni(t  -  2Ax)  )e"(bii/aii)2Ax 

+  4Ni(t  -  &T)e-(lDLi/*Li)&+   Ni(t)] 
where    At    =  data  time  increments 
and         A(t  -  2  At)  =  ZA(t)  evaluated  two  time  increments 

previously. 
The  value  z^(t)  is  evaluated  as 

Zi(t)  =  (1/aii)  (  At/3)  (ZAi(t)),  i  =  l,m. 
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To  obtain  the  values  of  the  system  response  for  a 
particular  time,  the  relationship 

(Y)  =  [*]  {z} 
where   {y}  =  system  response  vector  is  used,  using  the  column 
vector  [z]  for  the  time  at  which  results  are  desired.   The 
values  of  block  forces  caused  by  the  system  response  are 
calculated  every  0.02  seconds  for  the  first  30  seconds  of  the 
earthquake.   Only  the  first  3  0  seconds  are  needed  due  to  the 
nature  of  the  El  Centro  earthquake  the  seismic  acceleration 
after  3  0  seconds  are  small  when  compared  to  the  first  30 
seconds.   The  resultant  block  forces  and  vessel  motions  are 
then  checked  against  failure  criteria  established  in  Section 
3.   The  amplitude  of  the  seismic  acceleration  is  reduced  until 
no  block  failure  occurs.   Hence,  the  maximum  earthquake 
acceleration,  in  g's,  that  the  vessel-drydock  system  can 
withstand  without  failure  can  be  found. 

The  three  degree  of  freedom  vessel  response  computer 
program  is  listed  in  Appendix  1.   Since  all  three  linearized 
models  found  in  Equations  4.2.8,  4.3.11  and  4.4.11  are  in  the 
form 

[A]{y>  +  [B]{y}  =  (E(t)}, 
the  three  degree  of  freedom  computer  program  is  modified  to 
give  the  one  and  two  degree  model  results. 
7.1.2   Computer  Program  Input 

The  computer  program  listed  in  Appendix  1  evaluates  the 
horizontal  and  vertical  translations  and  rotation  of  the  three 
degree  of  freedom  vessel-drydock  system  in  increments  of  0.02 
seconds.   The  program  requires  the  following  input:   vessel 
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weight,  WEIGHT ,  in  kips,  keel  to  center  of  gravity  height 
(KG) ,  H,  in  inches,  vessel  mass  moment  of  inertia  about  the 
keel,  IK,  in  kip-in-sec2,  time  increments  of  ground 
acceleration  data,  DTAU,  in  seconds,  side  and  keel  piers 
vertical  and  horizontal  stiffnesses,  Ksv,  K^,  Kg^,  and  Kj^, 
in  kips/inch,  gravitational  constant,  GRAVITY,  in  inch/sec2, 
side  and  keel  pier  dimensions  of  width  and  height,  BASE  SIDE, 
BASE  KEEL,  HTSIDE  and  HTKEEL,  in  inches,  the  coefficient  of 
friction  between  the  block-block  and  ship-block  interfaces,  Ul 
and  U2,  horizontal  distance  between  the  center  of  vessel 
contact  with  port  and  starboard  side  piers,  BR,  in  inches, 
vertical/horizontal  ground  acceleration  ratio,  AMP, 
proportional  limits  of  side  and  keel  block  materials,  PLSIDE 
and  PLKEEL,  in  kips/in2,  total  cross-sectional  area  of  side 
and  keel  piers,  SIDE  AREA  and  KEEL  AREA,  in  in  2,  percent 
critical  damping,  ZETA,  vessel  hull  number,  HULL,  and  finally 
the  vessel-drydock  configuration  number,  NSYS.   The  following 
is  where  these  inputs  can  be  found: 
Table  5.1  NSYS  and  HULL 

Table  5.2  WEIGHT,  H,  Ik,  and  BR 

Table  5.3  BASESIDE,  BASEKEEL,  HTSIDE,  SIDEAREA 

and  KEELAREA 
Table  5.4  K^,  Ksv,  Kj^  and  Ksh. 

The  rest  of  the  inputs  are  assumed  to  be 
GRAVITY  =  384  in/sec2 
Ul  =  U2  =  0.5 
PLSIDE  =  PLKEEL  =0.7  kips/in2 
and     AMP  =0.5. 
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The  value  of  0.5  for  the  vertical/horizontal  ground  acceler- 
ation ratio  (AMP)  parameter  is  used  because  no  vertical 
component  of  the  El  Centro  earthquake  acceleration  was 
available  at  the  time  of  this  thesis. 

Damping  coefficients  equal  to  5%  of  critical  damping 
(i.e.  ZETA  =  0.05)  are  assumed  for  the  entire  system.   This 
analysis  for  damping  coefficients  can  be  found  in  Appendix  3 
and  is  implemented  directly  into  the  computer  program. 

Certain  idealized  system  failure  criteria  from  Section  3 
need  to  be  modified  to  accurately  represent  the  vessel-drydock 
system  as  it  actually  exists.   As  stated  in  Section  3.1,  block 
crushing  occurs  when  the  stress  on  the  block  exceeds  the 
block's  proportional  limit.   The  vessel  rests  upon  soft  wood 
caps  placed  on  top  of  the  keel  and  side  piers  when  in  the 
drydock  cradle.   Crushing  of  the  cap  is  not  considered  a 
failure  due  to  its  small  size.  For  the  blocking  system 
considered,  the  generated  stresses  are  transferred  through  the 
cap  to  the  top  drydock  block. 

Also,  a  side  pier  cannot  tip  inboard  because  the  vessel 
hull  is  not  physically  attached  to  the  pier  and  cannot  pull 
the  pier  beyond  it's  upright  position.   The  softwood  cap  on 
top  of  the  side  pier  causes  the  vessel  force  to  be  applied  12 
inches  from  the  inboard  edge.  As  shown  in  Figure  7.3,  this 
causes  the  static  force  vector  to  fall  outside  of  the  middle 
one-third  of  the  pier  base,  as  discussed  in  Section  3.   Since 
the  pier  will  only  tip  in  the  outboard  direction,  the 
resultant  vessel  force  must  fall  in  the  outboard  one-third  of 
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the  pier  base.   Thus,  this  failure  criteria  (from  Figure  7.3) 

for  the  side  pier  is 

HORIZONTAL  FORCE    2/3  Width  -  12.0   . 
VERTICAL  FORCE        Pier  Height 

7.1.3   Computer  Program  Testing 

In  order  to  validate  the  modal  anlaysis  method,  the  one 
degree  of  freedom  mathematical  model  is  checked  with  a  set  of 
sinusoidal  ground  accelerations  for  input.   The  input  ground 
motion  was 

E(t)  =  xg  =  323.95  sin(3.92t).  (7.1.1) 

This  sample  function  is  selected  because  its  magnitude  and 
frequency  closely  match  that  of  the  maximum  acceleration 
portion  of  the  El  Centro  earthquake.   Also,  vessel-drydock 
system  configuration  #1  is  chosen  for  this  initial  analysis. 
The  results  of  the  one  degree  of  freedom  case  will  be  compared 
to  the  closed  form  solution. 

The  one  degree  of  freedom  linear  equation  of  motion 

(Equation  4.2.7)  is  of  the  form 

',[  r       £      i  M  KG 

6  +  2  &>n  6  +  wn3e  = — —  E  (T  )  . 

k 
The  closed  form  solution  of  this  equation  is 

t 


e(t)   -"_*£  /  E(T)   e   ^      sina),(t-T)  di 

o 


•¥-££    [A(t)    sina),t    -    B(t)    cosco.t]  (7.1.2) 

I,  in),  a  a 

k   d 


where 


/ 


£  CO  T 

n 
(t)  eg-aj  ,  cosa)dT  di 

e 
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t  ?C0  T 

;  n 


and   B(t)  =  /  E(t)  ~-  -  sin  w^ 


o      e 

The  results  of  a  computer  run  with  vessel-drydock  system 
#1  and  sinusoidal  excitation  (Equation  7.1.1)  indicate  that 
the  maximum  rotation  angle  e  is  0.038243  radians  at  a  time  of 
2.02  seconds.   To  confirm  this  value,  Equation  7.1.2  has  to  be 
evaluated 


6(2.02)=-  VLEL       [A(2.02)    sino)d(2.02)    -B(2.02)    coswd(2  .  02)  ] 


2-02  Ku}    T 

/e 
323.95    sin(3.92x)  —  — — - —  cosco^    x    dx 
2 .  02^co  u 

0  e  n 

2.02 

/Cwnx 
323.95    sin(3.92x)    -77^ —  sinWdTdT* 

e  *     n 
0 


From  Section  6.3, 

wn  =  [((B2/2)KSV  -  W  KG)/IK)1/2 
=  5.580  rad/sec 
and   cod  =  wn(l  -  £2jl/2 

=  5.573  rad/sec. 
With  these  values  of  natural  and  damped  natural  frequencies, 
the  value  of  A (2. 02)  can  be  evaluated  using  the  relationship 


/  eat  sinbt   cosct  dt  =  eat  ^  sin (b-c) t- (b-c) cos (b-c) t 


2[a2+(b-c)23 


+   at  [a  sin  (b+c)  t-  (b+c)  cos  (b-t-c)  t] 
2[a2+(b+c)2] 
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Substituting  the  values 
t  =  2.02 


a  =  oon  =  .279 


b  =  3.92 
and   c  =  wd  =  5.573 
The  value  of  A (2. 02)  becomes 

323.95(e-5636)"1  *  -  .8307  =  -153.1629 
Similarly,  for  B(2.02), 
t 


at   .  ..        ,.     at  [ (b-c) sin (b-c) t+acos (b-c) t] 
e    sinbt  sinct  dt  =  e    -i-^— 


A  2[a2+(b-c)2] 

u  t 


at  [  (b+c) sin (b+c) t+acos (b+c) t] 


2[a2+(b+c)2] 
with  the  same  values  of  a,  b  and  c  as  before. 
Evaluating,  B(2.02)  becomes 

323.95  (e*5636)"1  *  -.2191  =  -40.3973 

Now, 

1  42    6979*193    0 

6(2.02)    =   -±-   *    *±z°212    -i^-u     [-153 . 1629*sin (2 . 02*5 . 573) 

2.54   2411000. *5. 573   -  (-40 . 3973*cos (2 . 02*5 . 573)  ] 

=  .038247  radians. 

Comparing  the  computer  generated  solution  (0.038243 
radians)  with  the  closed  form  solution  (.038247  radians),  it 
can  be  seen  that  the  computer  program  is  calculating  the  one 
degree  of  freedom  rotation  formula  correctly. 

One  other  simple  method  to  check  to  see  if  the  rotation 
displacement  time  history  is  correct  is  to  look  at  the 
waveform.   At  steady  state  condition,  6(t)  will  be  a 
sinusoidal  with  the  same  frequency  as  the  sinusoidal  ground 
acceleration.   An  amplitude  different  and  a  phase  shift  will 
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be  present.   The  steady  state  solution  [12]  of  a  equation  of 
motion  in  the  form 

mx  +   ex  +  kx  =  F  sinwt 


is 


x   = 


FJ 1 .     /        ^       -1    2£r    \ 

— L  =  -      sin  [cot-tan        ~  1 

k  V(l-r2)2+(2?r)2  \  l~r    / 


where  r  =  w/wn 

wn  =  (k/m) -1/2   =  Natural   Frequency 

K   =  1/2   c/  /  km  =  percent  critical  damping. 
In  this  case, 

wn  =  5.580   rad/sec 

r   =  3.92/5.580    =    .703 

5   =  .05 

|F|    =  M  KG — ^-=    1051012.9    lb. in 

1  '         2.54 

k  =  75071464  lb  in. 
so  now, 

6  steady  state  =  -027402  sin(3.92t  -  .1381)  radians  (7.1.3) 
Figure  7.1  represents  the  sinusoidal  ground  acceleration 
(Equation  7.1.1)  and  Figure  7.2  which  represents  the  rotation, 
6(t) ,  of  the  vessel-drydock  system  as  generated  by  the  modal 
analysis  computer  program.   Note  the  same  frequency  between 
Figures  7 . 1  and  7.2.   The  steady  state  portion  of  Figure  7 . 2 
is  matched  to  Equation  7.1.3.   In  conclusion  the  model 
analysis  method  correctly  predicts  the  one  degree  of  freedom 
mathematical  model. 
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7.1.4   Computer  Program  Results 

The  modal  analysis  method  can  now  be  used  with  the  El 
Centro  earthquake  data  to  predict  the  system  response.   Eleven 
vessel-drydock  configurations  (Section  5.0)  are  analyzed  for 
maximum  ground  acceleration  that  the  system  can  withstand 
without  block  failure.   Also,  the  first  block  failure  modes 
exhibited  by  the  systems  are  found. 

For  the  eleven  vessel-drydock  systems,  the  linearized 
one,  two  and  three  degree  of  motion  mathematical  model 
responses  are  calculated  and  listed  in  Table  7.3,  7.4,  and 
7.5,  respectively.   These  tables  can  be  found  on  pages  11 
thru 9-3.   The  program  further  calculates  that  vessel  liftoff 
(Section  3.4)  will  be  the  first  failure  mode  to  occur  when 
ground  accelerations  are  greater  than  the  maximum  permissible 
ground  acceleration  listed  in  Table  7.3,  7.4,  and  7.5.   The 
maximum  permissible  ground  acceleration  is  based  on  the 
reduction  of  the  El  Centro  acceleration  amplitude  until  no 
failure  occur.   However,  the  waveform  remains  unchanged. 

In  order  to  get  a  visual  picture  of  the  various 
responses,  vessel-drydock  system  #1  is  selected  to  be  plotted. 
Figures  7.4;  7.5a  and  b  and  7.6a,b,  and  c  are  the  appropriate 
responses  for  the  one,  two  and  three  degree  of  freedom  models. 
Only  the  first  thirty  seconds  of  the  earthquake  are  required 
to  give  the  maximum  system  response  as  seen  by  the  plots. 
Note  that  the  relative  horizontal  displacement  (Figure  7.6a) 
and  rotation  (7.6b)  time  histories  of  the  three  degree  of 
freedom  model  are  indentical  to  the  ones  found  in  the  two 
degree  of  freedom  model.   This  is  due  to  the  uncoupled 
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vertical  equation  of  motin  found  in  the  linearized  three 
degree  of  freedom  case. 

As  a  check  for  the  computer  generated  results,  the  system 
maximum  responses  due  to  full  magnitude  El  Centro  earthquake 
accelerations  are  compared  to  the  system  maximum  responses 
determined  by  response  spectrum  analysis  in  Table  7.1  and  7.2. 
These  tables  are  generated  by  combining  the  computer  results 
of  various  models  and  Tables  6.1  (Section  6.3)  and  A. 3. 3 
(Appendix  3) .   Since  the  response  spectrum  method  is  an 
approximate  one  (Section  6.3),  the  maximum  20%  error  between 
its  results  and  the  computer  results  in  the  one  degree  of 
freedom  case  is  reasonable.   Thus,  the  one  degree  program  is 
generating  accurate  results. 

In  the  two  degree  of  freedom  case,  the  response  spectrum 
method  predicts  the  worst  case  if  the  system's  two  natural 
modes  reach  a  maximum  at  the  same  time,  and  add  together. 
This  happened  in  the  rotational  case,  but  only  twice  in  the 
relative  displacement  case.   In  the  majority  of  the  systems, 
the  computer  generated  rotational  response  actually  exceeded 
the  rotational  response  predicted  by  modal  analysis.   The 
maximum  35%  difference  between  the  two  results  can  be 
accounted  for  because  of  the  approximate,  simplified  response 
spectrum  mentioned  earlier.   The  three  degree  of  freedom  case 
has  the  same  order  of  error  as  found  in  the  two  degree  case  as 
shown  in  Table  7.2.   Hence,  all  mathematical  models  perform 
satisfactory  using  the  modal  analysis  computer  program. 
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TABLE  7.1 

One  Degree  of  Freedom  Equation  of  Motion  Response  Comparison 
Modal  Analysis  Computer  Results  of  Linear  System 

I  e  | 

System  radians 

1  0.01474 

2  0.01618 

3  0.01436 

4  0.02966 

5  0.02966 

6  0.00832 

7  0.0115 

8  0.01712 

9  0.01512 
10  0.01714 
11 

Response  Spectrum  Method  Results  of  Linear  System 

I  e   I 
System  radians 

1  0.01594 

2  0.02619 

3  0.02149 

4  0.031 

5  0.031 

6  0.00698 

7  0.0068 

8  0.00987 

9  0.01445 

10  0.01299 

11  0.0145 
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TABLE  7.2 


Three  Degree  of  Freedom  Equation  of  Motion  Response  Comparison 
Modal  Analysis  Computer  Results  of  Linear  System 


System 


x 


(inches) 


1 

0.27971 

2 

0.22661 

3 

0.53288 

4 

0.58122 

5 

0.24721 

6 

0.36179 

7 

0.39118 

8 

0.39889 

9 

0.61065 

10 

0.57808 

11 

0.61730 

I   Y  I 
finches) 


07374 
10948 
15395 
18352 
11703 
06123 
0.06870 
0.07140 
0.11749 
0.11970 
0.11676 


( radians) 

0. 

,01458 

0. 

,01636 

0. 

,01546 

0. 

,03130 

0. 

,02992 

0. 

,00814 

0. 

,00879 

0. 

,01154 

0. 

,01587 

0. 

,01547 

0. 

,01734 

Response  Spectrum  Method  Results  of  Linear  System 


x 


System 


(inches) 


1 

0.33275 

2 

0.27141 

3 

0.65142 

4 

0.53788 

5 

0.21468 

6 

0.45271 

7 

0.46833 

8 

0.46597 

9 

0.64366 

10 

0.65864 

11 

0.67467 

1  y  1 

1  9    1 

finches) 

fradians) 

0.10859 

0.01646 

0.12225 

0.02348 

0.19089 

0.02350 

0.21491 

0.02350 

0.12980 

0.03202 

0.08025 

0.00683 

0.09708 

0.00866 

0.10762 

0.01009 

0.15381 

0.01552 

0.13662 

0.01342 

0.14955 

0.01552 

NOTE:   Two  Degree  of  Freedom  System  Maximum  Response 
correspond  to  x  and  Q . 
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7.2   Fourth  Order  Runge-Kutta  Analavsis  Solution 

7.2.1   Computer  Program  Development 
As  described  in  Section  6.2,  the  Runge-Kutta  scheme  uses 
a  set  of  first  order  differential  equations  to  describe  the 
second  order  vessel-drydock  system.   A  simple  example  of  this 
procedure  is  done  on  the  following  equation 

mx  +  ex  +  kx  =  f(t)  +  g(t)x 
let  r  =  x  and  rearranging  yields  two  first  order  equations, 
x  =  r 

r  =  -£r-£x  +  f(t)  +  g(t)x 
Now,  the  first  order  system  can  be  implemented  into  a  computer 
program. 

Briefly,  the  fourth  order  Runge-Kutta  numerical  scheme  is 
one  that  makes  use  of  predicted  velocity  and  acceleration  at  a 
given  time  to  compute  the  displacement  and  velocity  for  the 
next  time  increment.   This  can  be  done  by  the  basic 
relationships 

xn+l  =  xn  +  xn  At 
and       xn+l  =  *n  +  *n  At 

However,  x  and  x  are  not  found  by  differentiation  but 
through  the  manipulations  of  the  second  order  equation  into  a 
set  of  first  order  equations.   This  can  be  done  by  procedure 
shown  in  the  previous  paragraph. 

The  first  step  in  the  computer  program  development  is  to 
place  the  set  of  first  order  differential  equations  into  the 
Runge-Kutta  formulas  shown  in  Section  6.2.   The  evaluation  of 
these  formulas  will  predict  the  velocity  and  acceleration 
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that  will  be  used  in  the  next  time  increment.   Now,  the  system 
response,  both  displacements  and  associated  velocities,  can  be 
calculated  using  the  previous  information.   This  process  is 
marched  through  time  steps  in  our  to  produce  a  time  history  of 
the  system  response. 

Now,  the  values  of  block  forces  caused  by  the  system 
response  are  calculated  every  0.01  seconds  for  the  first  3  0 
seconds  of  the  earthquake.   The  first  3  0  seconds  are  only 
required  to  analyze  for  maximum  seismic  response.   The 
resultant  block  forces  and  vessel  motions  are  then  checked 
against  failure  criteria  found  in  Section  3.   The  amplitude  of 
the  seismic  acceleration  is  reduced  until  no  block  failure 
occurs.   Like  in  the  modal  analysis  method  (Section  7.1),  the 
maximum  earthquake  acceleration,  in  g's,  that  the 
vessel-drydock  system  can  withstand  without  failure  can  be 
found. 

The  Runge-Kutta  analysis  gives  the  ability  to  analyze 
both  the  linear  and  non-linear  equations  of  motion  for  the 
one,  two  and  three  degree  of  motion  mathematical  models.   The 
non-linear  three  degree  of  freedom  vessel  response  computer 
program  is  listed  in  Appendix  2.   The  program  is  developed 
around  Equations  4.4.5,  4.4.6  and  4.4.8  and  those  found  in 
Section  6.2.   Since  all  linear  and  non-linear  models  developed 
in  this  thesis  can  be  arranged  into  the  Runge-Kutta  formulas, 
the  non-linear  three  degree  of  freedom  comptuer  program  is 
modified  to  give  desired  model  results. 
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7.2.2  Computer  Program  Input 

The  computer  program  listed  in  Appendix  2  evaluates  the 
horizontal  and  vertical  translations  and  rotation  of  the 
non-linear  three  degree  of  freedom  vessel-drydock  system  in 
increments  of  0.01  seconds.   The  program  requires  the 
identical  input  found  in  Section  7.1.2  for  the  modal  analysis 
method.   Hence,  the  program  inputs  needs  not  be  listed  again. 
Like  in  the  previous  case,  Appendix  3  provides  the  damping 
coefficients  to  be  implemented  into  the  computer  program. 

7.2.3  Computer  Program  Testing 

In  order  to  validate  the  fourth  order  Runge-Kutta 
numerical  scheme,  the  linear  one  degree  of  freedom 
mathematical  model  is  checked  with  a  set  of  sinusoidal  ground 
accelerations  for  input.   The  same  input  ground  motion  that 
is  used  to  validate  the  modal  analysis  is  implemented  in  this 
case,  i.e. 

E(t)  =  xg  =  329.95  sin(3.92t).  (7.2.1) 

This  sample  function  is  selected  because  its  magnitude  and 
frequency  closely  match  that  of  the  maximum  acceleration 
portion  of  the  El  Centro  earthquake.   Also,  vessel-drydock 
system  configuration  #1  is  chosen  for  this  initial  analysis. 
The  results  of  the  linearized  one  degree  of  freedom  case  will 
be  compared  to  the  modal  analysis  and  closed  form  solutions. 

The  results  of  a  computer  run  with  vessel-drydock  system 
#1  and  sinusoidal  excitation  (Equation  7.2.1)  indicate  that 
the  maximum  rotation  angle  6  is  0.0382  3  6  radians  at  a  time  of 
2.02  seconds.   Comparing  this  the  computer  generated  solution 
with  the  modal  analysis  solution  (0.038243  radians)  and  the 
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closed  form  solution   (0.038247  radians)  found  in  Section 
7.1.3,  it  can  be  seen  that  the  computer  program  is  calculating 
the  linearized  one  degree  of  freedom  rotation  formula 
correctly.   Also,  Figure  7.7  which  represents  the  rotation,  e 
(t) ,  of  the  vessel-drydock  system  as  generated  by  the 
Runge-Kutta  computer  program  is  essentially  identical  to 
Figure  7.2  generated  by  the  modal  analysis  computer  program. 

In  conclusion,  the  fourth  order  Runge-Kutta  method 
correctly  predicts  the  linearized  one  degree  of  freedom 
mathematical  model.   Now,  the  computer  program  will  be  used 
for  the  linear  and  non-linear  one,  two  and  three  degree  of 
freedom  systems. 

7.2.4   Computer  Program  Results 

The  fourth  order  Runge-Kutta  analysis  can  be  implemented 
with  the  El  Centro  earthquake  data  to  predict  the  system 
response.   Eleven  vessel-drydock  configurations  (Section  5.0) 
are  analyzed  for  maximum  ground  acceleration  that  the  ystem 
can  withstand  without  block  failure.   Also,  the  first  block 
failure  modes  exhibited  by  the  systems  are  found. 

For  the  eleven  vessel-drydock  systems,  the  linear  and 
non-linear  one,  two  and  three  degree  of  motion  mathematical 
model  responses  are  calculated  and  listed  in  Table  7.3,  7.4 
and  7.5,  respectively.   These  tables  can  be  found  on  pages  ^1 
thru  9  3.   The  program  further  calculates  that  vessel  liftoff 
(Section  3.4)  will  be  the  first  failure  mode  to  occur  when 
ground  accelerations  are  greater  than  the  maximum  permissible 
ground  acceleration  listed  in  Table  7.3,  7.4  and  7.5. 
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In  order  to  get  a  visual  picture  of  the  various  linear 
and  non-linear  responses,  vessel-drydock  system  #1  is  selected 
to  be  plotted.   Figures  7.8,  7.9  a  and  b,  7.10  a,  b  and  c  are 
the  appropriate  response  for  the  linear  one,  two  and  three 
degree  of  freedom  models  and  Figures  7.11,  7.12a  and  b  and 
7.13  a,  b  and  c  are  for  the  non-linear  models.   Note  that 
these  figures  are  almost  identical  to  Figures  7.3,  7.4a  and  b 
and  7.5  a,  b  and  c  generated  by  the  modal  analysis  method 
(Section  7.1.5).   This  reinforces  the  equivalency  of  the 
results  found  by  the  two  numerical  shemes.   This  can  also  be 
seen  by  comparing  the  system  maximum  responses  of  the 
Runge-Kutta  method  to  that  of  the  modal  analysis  method.   Both 
are  listed  in  Table  7.3,  7.4  and  7.5.   The  maximum  error  for 
any  given  case  is  5%.   Hence,  all  mathematical  models  perform 
satisfactory  using  the  fourth  order  Runge-Kutta  computer 
program. 


-90- 


TABLE  7.3 
Modal  Analysis  Method  Results  of  Linear  System 


"max 

Time 

Maximum  Acceleration 

Svstem 

(radians) 

(sec) 

(cr's) 

1 

0.003390 

6.64 

0.0759 

2 

0.003883 

9.28 

0.0792 

3 

0.006033 

9.22 

0.1368 

4 

0.006822 

8.68 

0.0759 

5 

0.004513 

8.68 

0.0462 

6 

0.001999 

5.36 

0.0792 

7 

0.002422 

5.84 

0.0924 

8 

0.002760 

8.98 

0.0792 

9 

0.005138 

7.54 

0.0990 

10 

0.004537 

8.98 

0.0990 

11 

0.004972 

7.92 

0.0957 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Linear  System 


"max 

Time 

Maximum  Acceleration 

Svstem 

(radians) 

(sec) 

(cr's) 

1 

0.003390 

6.64 

0.0759 

2 

0.003882 

9.27 

0.0792 

3 

0.006032 

9.22 

0.1386 

4 

0.006822 

8.68 

0.0759 

5 

0.004153 

8.68 

0.0462 

6 

0.002004 

5.35 

0.0792 

7 

0.002426 

5.83 

0.0924 

8 

0.002763 

8.97 

0.0792 

9 

0.005141 

7.54 

0.0990 

10 

0.004543 

8.98 

0.0990 

11 

0.004972 

7.92 

0.0957 

4th  Order  Runge-Kutta  Numerical  Method  Results  of 

Non-Linear  System 


max 

Time 

Maximum  Acceleration 

Svstem 

(radians) 

(sec) 

fa's) 

1 

0.003391 

6.64 

0.0759 

2 

0.003883 

9.27 

0.0792 

3 

0.006032 

9.23 

0.1386 

4 

0.006825 

8.68 

0.0759 

5 

0.004154 

8.68 

0.0462 

6 

0.002004 

5.35 

0.0792 

7 

0.002427 

5.83 

0.0924 

8 

0.002764 

8.97 

0.0792 

9 

0.005139 

7.54 

0.0990 

10 

0.004546 

8.98 

0.0990 

11 

0.004970 

7.91 

0.0957 
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TABLE  7.4 
Two  Degree  of  Freedom  Equation  of  Motion  Response 
Modal  Analysis  Method  Results  of  Linear  System 


e 

max 

Maximum 

X 

max 

Time 

Time 

Acceleration 

System 

(inches)  ( 

seconds) 

(radians) 

(seconds) 

(e's) 

1 

0.064335 

8.38 

0.003355 

9.02 

0.0759 

2 

0.063873 

5.38 

0.003927 

9.26 

0.0792 

3 

0.207826 

5.42 

0.006033 

9.24 

0.1287 

4 

0.127868 

5.72 

0.006887 

8.66 

0.0726 

5 

0.032138 

5.80 

0.003891 

8.72 

0.0429 

6 

0.090448 

5.78 

0.002037 

5.36 

0.0825 

7 

0.109531 

5.46 

0.002463 

5.86 

0.0924 

8 

0.095734 

5.98 

0.002770 

7.82 

0.0759 

9 

0.195410 

6.06 

0.005080 

7.60 

0.1056 

10 

0.167644 

6.04 

0.004487 

7.84 

0.0957 

11 

0.172847 

6.08 

0.004856 

7.54 

0.0924 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Linear  System 


e 

max 

Maximum 

X 

max 

Time 

Time 

Acceleration 

System 

(inches)  ( 

seconds) 

(radians) 

(seconds) 

(p's) 

1 

0.064566 

8.38 

0.003358 

9.02 

0.0759 

2 

0.064934 

5.38 

0.003927 

9.26 

0.0792 

3 

0.207015 

5.56 

0.006034 

9.24 

0.1287 

4 

0.127317 

5.72 

0.006890 

8.66 

0.0726 

5 

0.032442 

5.80 

0.003891 

8.72 

0.0429 

6 

0.090984 

5.78 

0.002048 

5.36 

0.0825 

7 

0.110229 

5.46 

0.002466 

5.86 

0.0924 

8 

0.091996 

5.98 

0.002656 

7.81 

0.0759 

9 

0.194354 

6.06 

0.005074 

7.59 

0.1056 

10 

0.166760 

6.03 

0.004493 

7.83 

0.0957 

11 

0.174626 

6.07 

0.004860 

7.55 

0.0924 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Non-Linear  System 

Maximum 


max 

Time 

max 

Time  Ac 

celeration 

System 

(inches) 

(seconds) 

(radians) 

(seconds) 

(e's) 

1 

0.064597 

8.38 

0.003358 

9.02 

0.0759 

2 

0.064839 

5.38 

0.003926 

9.26 

0.0792 

3 

0.206443 

5.56 

0.006035 

9.24 

0.1287 

4 

0.127281 

5.72 

0.006892 

8.66 

0.0726 

5 

0.032442 

5.80 

0.003891 

8.72 

0.0429 

6 

0.090984 

5.78 

0.002048 

5.36 

0.0825 

7 

0.110274 

5.46 

0.002467 

5.86 

0.0924 

8 

0.092031 

5.98 

0.002656 

7.81 

0.0759 

9 

0.194361 

6.06 

0.005075 

7.59 

0.1056 

10 

0.166763 

6.03 

0.004493 

7.83 

0.0957 

11 

0.174636 

6.07 

0.004859 

7.55 

0.0924 
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TABLE  7.5 


Three  Degree  of  Freedom  Equation  of  Motion  Response 
Modal  Analysis  Method  Results  of  Linear  System 


e 

max 

Maximum 

max 

Time 

^max 

Time 

Time 

Acceleration 

System 

(inches) 

Csec .  ) 

Cinches) 

(sec . ) 

(radians) 

(sec .  ) 

(e's) 

1 

0.067132 

8.38 

0.017698 

5.68 

0.003501 

9.02 

0.0792 

2 

0.061211 

5.38 

0.025181 

8.06 

0.003765 

9.26 

0.0759 

3 

0.207826 

5.42 

0.060041 

5.64 

0.006033 

9.24 

0.1287 

4 

0.122056 

5.72 

0.038540 

5.64 

0.006574 

8.66 

0.0693 

5 

0.032138 

5.70 

0.015215 

5.38 

0.003891 

8.72 

0.0429 

6 

0.090448 

5.78 

0.019865 

5.32 

0.002037 

5.36 

0.0825 

7 

0.105619 

5.46 

0.018550 

6.68 

0.002375 

5.86 

0.0891 

8 

0.091745 

5.98 

0.016424 

5.36 

0.002654 

7.82 

0.0759 

9 

0.201516 

6.06 

0.038773 

5.80 

0.005238 

7.60 

0.1089 

10 

0.161863 

6.04 

0.033516 

5.38 

0.004332 

7.84 

0.0924 

11 

0.179020 

6.08 

0.033861 

5.40 

0.005030 

7.54 

0.0957 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Linear  System 


X 

max 

y 

max 

e 

Maximum 

Time 

Time 

max 

Time 

Acceleration 

System 

(inches) 

(sec .  ) 

(inches) 

(sec . ) 

(radians) 

(sec . ) 

(e's) 

1 

0.067373 

8.38 

0.017858 

5.37 

0.003504 

9.02 

0.0792 

2 

0.062228 

5.39 

0.025102 

8.06 

0.003761 

9.27 

0.0759 

3 

0.207015 

5.43 

0.059736 

5.64 

0.006304 

9.25 

0.1287 

4 

0.121529 

5.73 

0.038651 

5.65 

0.006577 

8.67 

0.0693 

5 

0.032442 

5.81 

0.014985 

5.39 

0.003891 

8.79 

0.0429 

6 

0.094623 

5.78 

0.019463 

5.77 

0.002130 

5.37 

0.0858 

7 

0.106292 

5.47 

0.018152 

6.69 

0.002378 

5.87 

0.0891 

8 

0.091996 

5.99 

0.016558 

5.37 

0.002656 

7.82 

0.0759 

9 

0.200428 

6.06 

0.039205 

5.40 

0.005233 

7.60 

0.1089 

10 

0.161009 

6.04 

0.033887 

5.39 

0.004338 

7.84 

0.0924 

11 

0.180863 

6.08 

0.034608 

5.40 

0.005033 

7.55 

0.0957 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Non- linear  System 


Maximum 

X 

max 

Time 

y 

-'max 

Time 

"max 

Time 

Acceleration 

System 

(inches) 

(sec . ) 

(inches) 

(sec . ) 

(radians) 

(sec . ) 

(p's) 

1 

0.067283 

8.38 

0.017858 

5.37 

0.003507 

9.02 

0.0792 

2 

0.062050 

5.39 

0.025102 

8.06 

0.003763 

9.27 

0.0759 

3 

0.206406 

5.57 

0.059736 

5.64 

0.006037 

9.25 

0.1287 

4 

0.121484 

5.73 

0.038651 

5.65 

0.006580 

8.67 

0.0693 

5 

0.032443 

5.81 

0.014985 

5.39 

0.003891 

8.79 

0.0429 

6 

0.094523 

5.78 

0.019463 

5.77 

0.002130 

5.37 

0.0858 

7 

0.106348 

5.47 

0.018152 

6.69 

0.002379 

5.87 

0.0891 

8 

0.091968 

5.99 

0.016558 

5.37 

0.002656 

7.82 

0.0759 

9 

0.200131 

6.06 

0.039205 

5.40 

0.005234 

7.61 

0.1089 

10 

0.160892 

6.04 

0.033887 

5.39 

0.004337 

7.84 

0.0924 

11 

0.180674 

6.08 

0.034608 

5.40 

0.005032 

7.55 

0.0957 
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7.3   Discussion  of  the  Models  and  Associated  Results 
Both  the  modal  analysis  and  the  Runge-Kutta  numerical 
methods  give  satisfactory  results  in  determining  the  system 
maximum  response  in  the  various  models  and  configurations. 
The  maximum  error  in  the  comparison  of  the  two  methods  for 
any  given  run  is  5%.   This  is  excellent  when  considering  that 
each  method  has  completely  different  basic  concepts  behind 
them.   The  modal  anlysis  method  runs  faster  to  obtain  the 
result  due  to  less  computer  operations.   But,  the  Runge-Kutta 
method  is  much  easier  to  program  into  the  computer  and  is  more 
flexible.   In  all,  both  methods  can  be  used  to  predict  the 
seismic  response  of  the  linearized  one,  two  and  three  degree 
of  freedom  models. 

In  the  comparison  of  the  one,  two  and  three  degree  of 
freedom  responses,  certain  conclusions  can  be  drawn.   The  one 
degree  of  freedom  model  need  only  be  used  if  system  maximum 
rotational  response  is  the  main  concern  to  drydock  system 
design.   This  model  gives  results  sufficiently  close  to  the 
rotation  generated  by  the  other  two  models.   Since  vessel 
liftoff  of  the  side  block  is  the  first  block  failure  to  be 
exhibited  by  the  vessel-drydock  system,  the  one  degree  of 
freedom  model  adequately  predicts  the  maximum  earthquake 
amplitude  the  system  can  withstand.   Most  of  the  vertical 
displacement  at  the  side  pier  is  due  to  rotation  and  not  of 
relative  vertical  displacement  as  found  in  the  three  degree 
model.   Hence,  the  one  degree  of  freedom  system  response 
should  be  used  to  determine  vessel  roation  and  maximum 
earthquake  amplitude. 
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The  three  degree  of  freedom  system  results  should  be  used 
for  a  parameter  study  on  block  failure  modes  and  for  block 
design.   This  is  because  only  the  three  degree  model  produces 
all  required  displacements  for  estimation  of  generated  forces 
in  both  keel  and  side  blocks.   However,  when  the  relative 
vertical  translations  generated  are  small,  the  two  degree 
model  can  be  used  in  block  design, .  The  one  degree  model  is 
truely  inadequate  for  block  design  since  crushing  and  vessel 
liftoff  of  the  side  blocks  are  the  only  failure  modes. 

The  non-linear  equations  of  motion  mathematical  models 
needs  to  be  addressed.  As  shown  in  Table  7.3,  7.4  and  7.5,  the 
non-linearity  of  the  system  caused  less  than  1%  difference 
between  the  linear  and  non-linear  responses.   Under  the  eleven 
respresentative  vessel-drydock  configurations,  the  non-linear 
terms  (ie.  y6  and  yg6)  should  be  removed  from  the  equations  of 
motion.   This  linearization  of  the  equations  of  motion  will  be 
verified  in  the  proceeding  discussion. 

The  next  step  in  the  studying  of  the  non-linearity  of  the 

system  is  to  determine,  when  the  non-linear  terms  mentioned 

earlier  should  be  included  in  the  equations  of  motion.   The 

non-linear  equation  of  motion  for  the  one  degree  of  freedom 

model  is 

c      k         —        — 
"     8  -    *1      -M  KG  ..     M  KG  ..  n 

6  +  IT  e  +  17  e   -IT"  xg  +  "T-  Yg  e 

k      k        k    ^      k    3 


Let  x„  =  a  sinwt 
YG  =  Y*g 
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M  KG  a/Ifc  =  2 

A2  =  kj/Ifc 

and  2  £A  =  Cg/I^. 

Substituting  and  rearranging  the  Equation  of  motion  gives 

6  +  2^9  +  A26  =  -8  sinwt(l  -y6)  .  (7.3.1) 

where  X  =  system  natural  frequency 

B  =  normalized  maximum  amplitude 
Y  =  ratio  of  vertical/horizontal  amplitude. 
Equation  7.3.1  not  only  represents  the  vessel-drydock  system 
but  any  other  system  than  can  be  reduced  to  this  form.   A 
parameter  study  is  performed  using  A,  B  and  Y  to  determine 
whether  or  not  to  use  non-linear  terms  in  the  equations  of 
motion.  The  selection  criterion  will  be  when  the  non-linear 
response  (ie.  Y  4   0)  is  different  to  the  linear  response  (ie  Y 
=  0)  by  10%.   If  for  a  given  X,   B  and  Y  this  criterion  is 
exceeded,  then  the  non-linear  term,Y6,  should  be  included.   If 
not,  exclude  them.   For  this  study,  let  co  =  3.92  rad/sec  and  t, 
=  0.05  to  emulate  the  sinusoidal  excitated  system  with  5% 
damping  discussed  in  Section  7.1.3.   Using  a  Runge-Kutta 
numerical  scheme,  Figures  7.14  a,  b,  c,  d  and  e  are  generated 
by  varying  the  three  parameters.   Note  that  a  negative  Y 
indicates  a  phase  shift  of  180°  between  the  vertical  and 
horizontal  sinusoidal  excitation.   The  linearity  selection 
criterion  of  10%  difference  holds  true  for  the  region  marked 
"N-L"  for  non-linear  for  a  given  B .   In  the  case  of  the  one 
degree  of  freedom  vessel-drydock  systems,  the  linearizing  of 
the  equations  of  motion  can  be  shown  on  Figure  7.14b  since  B  = 
0.5,Y=  0.5  and  A  ranges  from  2.870  rad/sec  to  10.254  rad/sec 
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(Table  6.1) . 

In  order  to  verify  the  linearity  selection  criterion,  a 
linear  point  selected  off  of  Figure  7.14b  of  3  =  0.5,  Y  =  5.58 
andY  =  2  0.0  is  implemented  into  the  Runge-Kutta  computer 
program  in  Section  7.2.   The  nonlinear  response  equals 
0.015839  radians  compared  to  the  linear  response  of  0.014741 
radians  when  using  the  El  Centro  earthquake  input.   This  is  a 
difference  of  7.5%  which  is  less  than  10%.   The  linearity 
selection  criterion  holds  true.   Now,  a  non-linear  point  of  3 
=0.5,  X  =  8.327  and  y  =  15.0  is  implemented  into  the 
Runge-Kutta  scheme.   The  non-linear  and  linear  responses  are 
.009652  and  .008664,  respectively.   This  is  a  difference  of 
11.4%  which  is  greater  than  10%.   The  linearity  selection 
criterion  is  validated.   Thus,  vessel-drydock  system  should  be 
checked  before  selection  of  linear  or  non-linear  models. 

As  mentioned  earlier,  Equation  7.3.1  can  be  used  to 
validate  systems  other  than  the  vessel-drydock  systems 
subjected  to  sinusoidal  ground  acceleration  whose  frequency 
matches  that  of  the  El  Centro  earthquake,  ie.  w  =  3.92 
rad/sec.   this  can  be  done  by  modifying  the  appropriate 
linearity  selection  criteria  figure  with 

\l    =  A/3.92 
and  replotting  the  figure.   This  essentially  changes  Equation 
7.3.1  into  a  non-dimensional  of 

0-L  +  2^X16!  +  Xx2  6-l  =  -&!  sinti  (1  -  Yi©!) 

where         \  =    Q/^2 

•         • 

e1  =  e/oj 
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1 

xl 

= 

X  /  03 

61 

= 

6/co2 

\ 

= 

Y 

and  in  this  case,  w  =  3.92  rad/sec. 
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FIGURE  7.1 
Sinusoidal  Ground  Acceleration 
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FIGURE  7.2 
Vessel  Seismic  Response  in  Rotation  using  Modal  Analysis  with 
Sinusoidal  Ground  Acceleration 
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FIGURE  7.3 
Side  Pier  Forces,  with  Softwood  Cap  Block 
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FIGURE  7.4 
Vessel  Seismic  Response  in  Rotation  for  One  Degree  of  Freedom 

Model  using  Modal  Analysis 
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FIGURE  7.5a 
Vessel  Seismic  Response  in  Relative  Horizontal  Translation  for 
Two  Degree  of  Freedom  Model  using  Modal  Analysis 
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FIGURE  7.5b 
Vessel  Seismic  Response  in  Rotation  for  Two  Degree  of  Freedom 

Model  using  Modal  Analysis 
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FIGURE  7.6a 
Vessel  Seismic  Response  in  Relative  Horizontal  Translation  for 
Three  Degree  of  Freedom  Model  using  Modal  Analysis 
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FIGURE  7.6b 
Vessel  Seismic  Response  in  Relative  Vertical  Translation  for 
Three  Degree  of  Freedom  Model  using  Modal  Analysis 
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FIGURE  7.6c 
Vessel  Seismic  Response  in  Rotation  for  Three  Degree  of 
Freedom  Model  using  Modal  Analysis 
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FIGURE  7.7 
Vessel  Seismic  Response  in  Rotation  using  Runge-Kutta 
Numerical  Analysis  with  Sinusoidal  Ground  Acceleration 
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FIGURE  7.8 
Vessel  Seismic  Response  in  Rotation  for  One  Degree  of  Freedom 
Model  using  Runge-Kutta  Analysis 
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FIGURE  7.9a 
Vessel  Seismic  Response  in  Relative  Horizontal  Translation  for 
Two  Degree  of  Freedom  using  Runge-Kutta  Analysis 
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FIGURE  7.9b 
Vessel  Seismic  Response  in  Rotation  for  Two-Degree  of  Freedom 
Model  using  Runge-Kutta  Analysis 
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FIGURE  7.10a 
Vessel  Seismic  Response  in  Relative  Horizontal  Translation  for 
Three  Degree  of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.10b 
Vessel  Seismic  Response  in  Relative  Vertical  Translation  for 
Three  Degree  of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.10c 
Vessel  Seismic  Response  in  Rotation  for  Three  Degree  of 
Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.11 
Vessel  Seismic  Response  in  Rotation  for  Non-Linear  One  Degree 
of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.12a 
Vessel  Seismic  Response  in  Relative  Horizontal  Translation 
for  Non-linear  Two  Degree  of  Freedom  Model  using  Runge-Kutta 

Analysis 
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FIGURE  7.12b 
Vessel  Seismic  Response  in  Rotation  for  Non-linear  Two  Degree 
of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.13a 
Vessel  Seismic  Response  in  Relative  Horizontal  Translation 
for  Non-linear  Three  Degree  of  Freedom  Model  using  Runge- 

Kutta  Analysis 


-118- 


(/J 

Ld 
I 

Z20E-3 


u 
z 
u 
o 
< 
-J 

Cl 
V) 

I — I 

Q 


< 


DC 
Ld 

> 

Ld 
> 


1 


E= 


10 

5 


-  10 

-  15 


0    ^^--J|| 

5   - 


^20E-3 
< 

-J  0 

u 

DC 


5 


10  15  20 

TIME    IN   SECONDS 


25 


30 


FIGURE   7.13b 

Vessel  Seismic  Response  in  Relative  Vertical  Translation  for 

Non-linear  Three  Degree  of  Freedom  Model  using  Runge-Kutta 

Analysis 
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FIGURE  7.13c 
Vessel  Seismic  Response  in  Rotation  for  Non-linear  Three 
Degree  of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.14b 
Linearity  Selection  Criteria,  5%  Damping,  Beta 
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Linearity  Selection  Criteria,    5%   Damping,    Beta  =1.0 
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Linearity  Selection  Criteria,  5%  Damping,  Beta 
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FIGURE  7.l4e 
Linearity  Selection  Criteria,  5%  Damping,  Beta  =10.0 
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8.0   SUMMARY.  CONCLUSIONS.  AND  RECOMMENDATIONS  FOR  FURTHER 

STUDY 

8. 1  Summary 

This  thesis  investigates  four  models  of  the  vessel- 
drydock  system  for  seismic  analysis.   These  models  are 
implemented  to  determine  vessel  response  to  seismic  loading 
for  eleven  representative  submarine-drydock  systems.   Using 
Equation  4.1.4,  the  quasi-static  force  model  generates  various 
vessel  responses  listed  on  Tabe  6.2.   For  example,  the 
submarine  in  system  #1  will  rotate  about  its  keel  0.00809  when 
subjected  to  a  strong  earthquake.   The  one  degree  of  freedom 
model  predicts  a  maximum  vessel  rotation  for  system  #1  of 
0.01471  radians  in  response  to  the  El  Centro  earthquake  ground 
acceleration  history.   This  rotation  of  0.01471  radians  is 
verified  by  the  two  and  three  degree  of  freedom  models.   See 
Tables  7.2  and  7.3.   This  maximum  seismic  response  of  system 
#1  can  also  be  confirmed  by  the  response  spectrum  method  as 
described  in  Section  6.3  and  reference  [1].   The  response 
spectrum  method  gives  the  rotation  as  0.01594  radians.   All 
eleven  systems  have  similar  differences  between  the  quasi- 
static  force  model  and  the  various  degree  of  freedom  models. 

More  important  than  the  maximum  vessel  response  due  to 
the  El  Centro  earthquake  of  .33g  magnitude  is  the  maximum 
permissible  earthquake  acceleration  that  a  vessel-drydock  can 
withstand  without  failure.   The  maximum  permissible 
accelerations  and  associated  vessel  responses  for  the  one,  two 
and  three  degree  of  freedom  are  listed  in  Table  7.3,  7.4  and 
7.5,  respectively.   For  example,  the  one  degree  of  freedom 
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model  for  system  #1  predicts  a  maximum  permissible 
acceleration  of  0.0759  g's  with  maximum  vessel  rotation  of 
0.00339  radians.   The  two  degree  of  freedom  determines  a 
maximum  permissible  acceleration  equalling  .0759  g's  with 
maximum  vessel  rotation  and  relative  horizontal  displacement 
of  0.003358  radians  and  0.064566  inches,  respectively. 
Finally,  the  maximum  permissible  acceleration  of  0.0792  g's 
with  associated  vessel  responses  (rotation,  relative 
horizontal  and  vertical  displacements)  of  0.003504  radians, 
0.067373  inches,  and  0.017858  inches.   The  slightly  higher 
maximum  permissible  acceleration  found  in  three  degree 
model  is  due  to  the  beneficial  effect  of  the  vertical 
displacement  delaying  failure  of  the  system.   However,  as 
shown  by  the  other  systems,  this  is  not  always  the  case.   In 
summary,  the  range  of  magnitudes  known  as  maximum  permissible 
accelerations  that  the  eleven  systems  could  withstand  without 
liftoff  is  13%  to  42%  of  the  magnitude  of  the  El  Centro 
earthquake  depending  on  system  and  model  used. 

8 . 2   Conclusions 

The  quasi-static  method  currently  used  by  the  U.S.  Navy 
for  seismic  response  analysis  underestimates  the  block  forces 
caused  by  an  earthquake  the  magnitude  of  the  El  Centro 
earthquake.   The  one  degree  of  freedom  model  adequately 
predicts  the  vessel  rotational  response  to  seismic  loading  but 
is  ineffective  in  the  analysis  of  possible  failures  to  the 
vessel-drydock  system.   The  two  degree  of  freedom  model  does  a 
good  job  of  predicting  vessel  rotational  and  horizontal 
translation  response  to  seismic  loading  provided  the  response 
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due  to  vertical  acceleration  is  small,  which  is  the  case  for 
the  eleven  systems  analyzed.   Finally,  the  three  degree  of 
freedom  model  predicts  the  vessel  rotational,  horizontal,  and 
vertical  translation  responses  and  can  be  used  to  completely 
investigate  various  types  of  failure  to  the  vessel-drydock 
system  in  all  eleven  vessel-drydock  systems  examined.   The 
vessel  would  experience  side  pier  liftoff  failure  during  an 
earthquake  with  the  magnitude  of  the  El  Centro  earthquake  (ie. 
0.33  g's).   The  various  vessel-drydock  systems  would  not 
remain  intact  during  this  magnitude  of  earthquake.   Also  noted 
that  the  fourth  order  Runge-Kutta  numerical  scheme  should  be 
used  to  evaluate  system  response  because  of  its  simplicity, 
results  accurancy,  and  ease  of  introducing  terms  into  the 
Equation  of  motion. 

8. 3   Recommendations  for  Further  Study 

The  seismic  response  of  drydocked  vessels  needs  further 
investigation  in  three  areas.   First,  a  three  degree  of 
freedom  model  which  will  allow  vessel  liftoff  should  be 
studied  so  that  other  failure  mode  will  occur.   Second,  a 
study  of  the  vessel-drydock  system  parameters  needs  to  be 
carried  out  to  explore  ways  in  which  drydocked  vessel  seismic 
response  can  be  decreased.   Third,  an  in  depth  study  of  the 
drydock  block  itself  in  order  to  provide  failure  modes  along 
with  stiffness  and  damping  characteristics  for  implementation 
into  mathematical  models. 
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APPENDIX  1 


C       LINEARIZED  THREE  DEGREE  OE  FREEDOM  SYSTEM  RESPONSE  USING 
C        MODAL  ANALYSIS  METHOD 

integer  ia,ib,lc,l,mm,n,i;Job,  iz,ier  ,hull,nsys,  flaglO 
integer  flagl, flag2, flag3, flag4, flag5, flag6, flag7, flag8 
real  a (6,6) ,b (6,6) ,beta,wk(72) , weight, h, Ik, gravity 
real  m(3,3) ,cx(3,3)  k(3,3) , sldearea,keelarea,plside 
real  bases lde,basekeel ,htside,htkeel ,critl,crit2,crit3 
real  an  (6,6)  ,  ac  (6000)  , dtau,maxx,maxt,maxy , timex,  timet 
real  rfl,r f2,rf3,hfl,hf2,hf3,ampacc,mass,ampacmax 
real  kvs , kvk , khs , khk , base , ht , counter 

real  timel , tlme2 , time 3, tlme4, time5, time6 , time7 , time8 
real  x(6000)  ,t  (6000)  ,y(6000)  ,ri(7)  ,sl(7)  ,pi(7)  ,XSCL(6) 
real  bbb,ccc, wl2, wl , w22 , w2 ,w32,w3, model, mode 3, capwidth 
real  mmxl ,mmangl ,mmx3,mmang3,br ,amp,plkeel,crit4,ul,u2 
real  timey ,mmmmml ,mmmmm2 ,mmmmm3,mmmmm4 
complex  alfa(6)  ,z(6,6)  , ad  (6,6)  ,bd(6,6)  ,aa(6,6)  ,bb(6,6) 
complex  g(6) ,v(6) ,yy (6,6000) ,ABC(6) ,zt(6,6) ,betal 

CHARACTER* 40  XLABEL, YLABEL, YYLABEL, YYYLABEL,DEC 


C  READ  IN  VESSEL  AND  DRYDOCK  DATA;  VESSEL  WEIGHT, KG, I (ABOUT  KEEL) 

C  TIME  INCREMENT  OF  DATA  POINTS, VERTICAL  STIFFNESS  OF  SIDE  AND 

C  KEEL  PIERS, HORIZONTAL  STIFFNESS  OF  SIDE  AND  KEEL  PIERS, 

C  GAVITATIONAL  CONSTANT, SIDE  BLOCK  BASE  AND  HEIGHT, 

C  KEEL  BLOCK  BASE  AND  HEIGHT, 

C  BLOCK-BLOCK  AND  BLOCK-HULL  FRICTION  COEFFICIENTS, 

C  SIDE  AND  KEEL  BLOCK'S  PROPORTIONAL  LIMIT, 

C  SIDE  PIER-VESSEL  CONTACT  AREA, KEEL  PIER-VESSEL  CONTACT  AREA,   - 

C  CAP  BLOCK  INCLINATION  ANGLE. 

read (44, *)  weight, h, Ik, dtau, kvs, kvk, khs, khk, gravity 
read(44,*)  baseslde,  basekeel ,htside,htkeel,ul,u2 
read  (44, *)  br , amp,plside,plkeel ,sidearea,keelarea,zeta 
read (44,*)  hull,nsys 

write  (6,*)  'do  you  want  response  plots?  (y  or  n) ' 
read  (5, 15)  dec 
15      format  (a) 
do  10,1=1,6 
do  ll,j=l,6 
a(l,  j)=0.0 
b(i,  j)=0.0 

11  continue 
10      continue 

do  12,1=1,3 
do  13,J=1,3 
m(l,  j)=0.0 
k(i,  j)=0.0 
ex  (1,3)  =0.0 
13      continue 

12  continue 
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CALCULATE  SYSTEM  PARAMETERS 

mass=weight/gravity 

beta=asln(sqrt(br**2/(4*h**2) ) ) 

m(l,l)  =mass 
m(l ,  3) =h*mass 
m(2,  2)=mass 
m (3, 1) =mass*h 
m(3,3)=Ik 

k(l,l)  =  (2*khs+khk) 

k(2,2)=(2*kvs+kvk) 

k(3,3)  =  ((0.5*kvs* (br**2 .0) ) - (weight *h) ) 

DETERMINE  NATURAL  FREQUENCIES  OF  SYSTEM 

bbb=-(m(l,l)  *k(3,3)+m(3,3)  *k(l,l)-m(l,3)*k(3,l)-m(3,l)  *k(l,3)) 

/(m(l,l)*m(3,3)-m(l,3)*m(3,l)) 

ccc=(k(l,l)  *k(3,3)-k(l,3)*k(3,l))/(m(l,l)*m(3,3)-m(l,3)*m(3,l)) 

NATURAL  FREQ.    MODE   #1 

wl2=(-bbb-sqrt(bbb**2.0-4*ccc))/2 
wl=sqrt (wl2) 

NATURAL  FREQ.  MODE  #2 

v22=k(2,2)/m(2,2) 

w2=sqrt (w22) 

NATURAL  FREQ.  MODE  #3 

w32=(-bbb+sqrt(bbb**2.0-4*ccc) ) /2 
w3=sqrt (w32) 

MODE  SHAPE  #1  &  #3 

model=(m(l,3)  *wl2-k  (1,  3) )  /(-m(l,  1)  *wl2+k(l,l)) 

mode3=(m(l,3)*w32-k(l,3))/(-m(l,l)*w32+k(l,l)) 

DETERMINE  C11,C13,C31,C33 

mmxl=m (1, 1) +m(l, 3) /model 

mmangl=model*m (3, 1)  +m(3,  3) 

mmv^rmi  (l,  l)  +m  (1 ,  3)  /mode3 

mmang3=mode3*m (3, 1) +m (3, 3) 

mmmmml=2*zeta*mmxl*wl 

mmmmm2=2 * zeta *mmx 3 * w 3 

mmmmm3=2  *  zeta  *mmangl  *  wl 

mmmmm4=2*zeta*mmang3*w3 


cx(l, 3) = (mmmmml-mmmmm2) / (l/model-l/mode3) 
ex  (1, 1) =mmmmml- (ex (1, 3) /model) 
cx(2,2)=2*zeta*m(2,2) *w2 
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ex  (3, 1)  =  (mmmmm3-mmmmin4)  /  (model-mode3) 
ex  (3,  3)  =mmmmm3-  (ex  (3,1)  *model) 

SET  UP  A  AND  B  MATRICES, PERFORM  EIGENVALUE  PROBLEM,  AND 
DECOUPLE  EQUATIONS  OF  MOTION 


a(l,4)  = 

m(l,l) 

a(l,6)  = 

m(l,3) 

a(2,5)  = 

m(2,2) 

a(3,4)  = 

m(3,l) 

a(3,6)  = 

m(3,3) 

a(4,l)  = 

0(1,1) 

a  (4,  3)  = 

m(l,3) 

a(4,4)  = 

cx(l,l) 

a(4,6)  = 

ex  (1,3) 

a(5,2)  = 

m(2,2) 

a  (5,5)  = 

cx(2,2) 

a  (6,1)  = 

m(3,l) 

a(6,3)  = 

m(3,3) 

a(6,4)  = 

cx(3,l) 

a(6,6)  = 

cx(3,3) 

b(l,l)  = 

-m-(l.l) 

b(l,3)  = 

-m(l,3) 

b(2,2)  = 

-m(2,2) 

b(3,l)  = 

-m(3,l) 

b(3,3)  = 

-m(3,3) 

b(4,4)  = 

k(l,l) 

b(4,6)  = 

k(l,3) 

b(5,5)  = 

k(2,2) 

b(6,4)  = 

k(3,l) 

b(6,6)  = 

k(3,3) 

do  100, 

i=l,6 

do  110, 

3=1.6 

an(i,j)  = 

=  -1.0*a(i,j) 

110 

continue 

100 

continue 

201 
200 


ia=6 

ib=6 

ic=6 

iz=6 

n=6 

ljob=2 

call  elgzf (b, ia,an, lb,n, i job,alfa,betal,z, iz,wk, ier) 

WRITE  (6,*)  wk(l)  ,ler 

DO  200,1=1,6 
do  201»j=l,6 

Zt(j,l)=  2(1, J) 

continue 

CONTINUE 

call  mult  (zt,a,aa) 

call  multc  (aa,z,ad) 

call  mult  (zt,b,bb) 

call  multc  (bb,z,bd) 

do  204,1=1,6 

v(i)=-1.0*bd(i,l)/ad(l,l) 
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204 


continue 


C       READ  IN  ACCELERATION  DATA 

do  300,n=l,5001,5 

read  (45,*)  ac  (n) , ac  (n+1) ,ac  (n+2) ,ac (n+3) ,ac (n+4) 
300     continue 
C       ESTABLISH  FAILURE  CRITERIA  AND  FLAGS 


10000 


50000 


crltl=  mln (ul, (u2*cos (beta) +sin (beta) ) / 

(cos (beta) -u2*sin (beta) ) ) 

crit2=min  (ul,u2) 

crit3=  (0 .66*bassside-12 .0) /htside 

crit4=basekeel/(6*htkeel) 

ampacc=l .0 

counter=0 .0 

ampacmax=0 .0 

continue 

flagl=0 

flag2=0 

flag3=0 

flag4=0 

flag5=0 

flag6=0 

flag7=0 

flag8=0 

flagl0=0 

do  50000  1=1,6 

ABC(i)=0.0 

continue 

maxx=0 . 0 

maxt=0.0 

maxy=0 . 0 

mm=0 

SOLVE  FOR  Y,THE  COLUMN  MATRIX  WHICH  IS  THE  SOLUTION  IN 

THE  IMAGINARY  COORDINATE  SYSTEM  WHERE  THE  EQUATIONS  OF  MOTION 

ARE  DECOUPLED 


302 


do  301,1=3,3501,2 

do  302,1=1,6 

g(i)  =  (ABC(i)+ac(l-2)/2.54)  *exp(v(l)  *2*dtau) 

+4.0*  (ac(l-l)/2.54)  *exp(v(i)  *dtau)  +ac  (1)  /2  .54 

ABC(l)=g(l) 

yy(l,l)  =  ((zt(l,4)+amp*zt(i,5)+zt(i,6)  *h)  *mass/ad (1, i) )  *g(i) 

*  ampacc  * -dtau/3 

continue 

mm=mm+l 


C       USING  THE  MODAL  MATRIX  (  4th, 5th,  AND  6th  ROWS  ),  OBTAIN 
C        VALUES  FOR  TRANSLATIONS  AND  ROTATION 

do  303,  1=1,6 
si(i)=z(4,i)*yy(i,l)+si(i-l) 

303     continue 

x(nrai)  =si  (6) 

if  (abs (x (mm) ) .gt .abs (maxx) )  then 

timex=dtau* (1-1) 

maxx=x (mm) 
endif 
do  304,  1=1,6 
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ri(i)=z(6,i)*yy(i,l)+rl(i-l) 

304  continue 

t  (mm)=rl  (6) 

if  (abs (t (mm) ) .  gt .abs (maxt) )  then 

tlmet=dtau* (1-1) 

maxt=t (mm) 
endlf 

do  305,  1=1,6 
pi(i)=2(5.i)*yy(i,l)+pl(i-l) 

305  continue 
y(mm)=pi(6) 

if  (abs (y (mm) ) .gt .abs (maxy) )  then 

timey=dtau* (1-1) 

maxy=y (mm) 
endif 

C       CALCULATE  VERTICAL  AND  HORIZONTAL  EORCES  CAUSED  BY  VESSEL, 
C        TEST  EOR  FAILURE 

C       CALCULATE  FORCES  ON  SIDE/KEEL  BLOCKS 

rfl=kvs*  (  (veightA(2,2))-y(mm)-(br/2)  *t(mm)) 

rf2=kvs*  (  (welght/k(2,2)  ) -y  (mm)  +  (br/2)  *t(mm)) 

rf3=kvk*  ((weightA(2,2))-y(mm)) 

hfl=khs* (x(mm)) 

hf2=khs* (x(mm)) 

hf3=khk*x(mm) 
C       TEST  FOR  SIDE  BLOCK  SLIDING 

If  (flagl.eq.l)  then 

go  to  400 
else  if  (hfl.lt.O.O.and.rfl.gt.O.O.and.abs (hfl/rfl) .gt.crltl)  then 

timel=  dtau*  (1-1) 

flagl=l 
else  if  (hf2.gt.0.0.and.rf2.gt.0.0.and.abs (hf2/rf2) .gt.critl)  then 

timel=dtau*  (1-1) 

flagl=l 
endif 
xl=x (mm) 
yl=y (mm) 
tl=t (mm) 
400     continue 

C       TEST  FOR  KEEL  BLOCK  SLIDING 

if  (flag2.eq.l)  then 

go  to  410 
else  if  (rf3.gt.0.0.and.abs (hf3/rf3) .gt.crit2)  then 
time2=dtau* (1-1) 
flag2=l 
endlf 
x2=x (mm) 
y2=y (mm) 
t2=t (mm) 
410     continue 
C       TEST  FOR  SIDE  BLOCK  OVERTURNING 

if  (flag3.eq.l)  then 

go  to  420 
else  If  (hfl.lt. CO. and. rfl.gt. 0.0. and. abs(hfl/rfl)  .gt.crlt3)  then 

time3=  dtau* (1-1) 
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flag3=l 
else  if  (hf2.gt.O.O.and.rf2.gt.O.O.and.abs (hf2/rf2) .gt.crlt3)  then 

time3=dtau* (1-1) 

flag3=l 
endlf 
x3=x(mm) 
y3=y  (mm) 
t3=t (mm) 
420     continue 

C       TEST  FOR  KEEL  BLOCK  OVERTURNING 

if  (flag4.eq.l)  then 

go  to  430 
else  if  (rf3.gt.0.0.and.abs (hf3/rf3) .gt.crlt4)  then 

time4=dtau* (1-1) 

flag4=l 
endif 
x4=x (mm) 
y4=y (mm) 
t4=t (mm) 
430     continue 

C       TEST  FOR  SIDE  BLOCK  LIFTOFF 

if  (flagS.eq.l)  then 

go  to  440 
else  if  (rfl.lt. 0^0  .or.  rf2.lt. 0.0)  then 

time5=dtau*  (1-1)  "• 

flag5=l 
endif 
x5=x (mm) 
y5=y (mm) 
t5=t (mm) 
440     continue 

C       TEST  FOR  KEEL  BLOCK  LIFTOFF 

if  (flag6.eq.l)  then 

go  to  450 
else  if  (rf3.lt. 0.0)  then 

time6=dtau*  (1-1) 

flag6=l 
endif 
x6=x (mm) 
y6=y (mm) 
t6=t (mm) 
450     continue 

C       TEST  FOR  SIDE  BLOCK  CRUSHING 

if  (flag7.eq.l)  then 

go  to  460 
else  if  (rfl.gt.0.0  .and.  (rfl/sldearea) .gt .plside)  then 

flag7=l 

tlme7=dtau* (1-1) 

else  if  (rf2.gt.0.0  .and.  (rf2/sidearea) .gt. plside)  then 
flag7=l 
time7=dtau*  (1-1) 
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460 
C 


endif 
x7=x (mm) 
y7=y (mm) 
t7=t (mm) 

continue 

TEST  FOR  KEEL  BLOCK  CRUSHING 


470 
301 


if  (flag8.eq.l)  then 

go  to  470 
else  if  (rf3.gt.0.0  .and.  (r f 3/keelarea) .gt .plkeel)  then 

flag8=l 

time8=dtau* (1-1) 
endif 
x8=x (mm) 
y8=y (mm) 
t8=t (mm) 
continue 

continue 


C       PLOT  RESULTS 

go  to  999 
60000   continue 

if (DEC.EQ. 'N')  THEN 

write  (6,*)  'I  am  finishing.' 
GO  TO  998 

endif 

write  (6,*)  'I  am  plotting.' 

XSCL (1) =0.0 

XSCL (2) =30.0     -&■ 

XLABEL='TIME  IN  SECONDS' 

YLABEL= 'ROTATION  IN  RADIANS' 

YYLABEL= 'RELATIVE  HORIZONTIAL  DISPLACEMENT  IN  INCHES* 

YYYLABEL= 'RELATIVE  VERTICAL  DISPLACEMENT  IN  INCHES' 

CALL  QPICTR  (X, 1 , 1500, QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) 
+   QYLAB (YYLABEL) , QLABEL (4)) 

CALL  QPICTR  (T, 1 , 1500, QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) 
+   QYLAB (YLABEL) , QLABEL (4) ) 

CALL  QPICTR  (Y, 1, 1500, QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) 
+   QYLAB (YYYLABEL) , QLABEL (4)) 

998  go  to  20000 

999  CONTINUE 


4000 
4050 
4100 
4150 
4200 

4250 


12, lx, *****') 


if (ampacc.eq. 1 .0)  then 

write (46 ,4000)  nsys 
format (lx,/,28x, ' ****  System 
write (46,4050)  hull 

format (lx,/, 30x, '**  Hull  ' , 13, lx, ' ** ' ) 
write (46,4100) 

format (lx,//,28x, ' *  Ship  Parameters  *') 
write  (46,4150) 

format (lx,/,5x, 'Weight' ,8x, 'Moment  of  Inertia' ,9x, 'K.G. ' ) 
write  (46, 4200)  weight, Ik, h 

format (lx, f 9.1, lx, 'kips' , lx, fll .1, lx, 'klps-in-sec2 ' , 
+3x,f6.1,lX,'ins') 
write (46, 4250) 
format (lx,//,26x, '*  Drydock  Parameter  *') 
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write (46 ,4300) 
4300    format (lx,/,lx, 'Side  Block  Height ', 3x, 'Side  Block  Width", 
+3x,'Keel  Block  Height ', 3x, 'Keel  Block  Width') 

write (46,4350)  htslde, bases lde,htkeel ,basekeel 
4350    format (2x,f 6. l,lx, ' ins ' , llx, f6. 1, lx, ' ins ' , llx, f6. 1, lx, 'ins' , 
+9x,f6.1,lx, 'ins') 

write  (46,4400) 
4400    format (lx,/, lx, ' Side-to-Side  Pier  Distance*) 

write (46, 4450)  br 
4450    format (Ix,t7,f6.1,lx, 'ins') 

write (46, 4470) 
4470    format (lx,/, '  Total  Side  Pier  Contact  Area' 
+,3x, 'Total  Keel  Pier  Contact  Area') 

write (46,4475)  sidearea,keelarea 
4475    format (lx,8x,f 11. l,lx, 'in2' , 14x, fll . 1 , lx, 'ln2') 

write (46,4500) 
4500     format (lx,/,lx, 'Block-on-Block  Friction  Coef f ' , 3x, 'Hull-on-Block 
+  Friction  Coeff) 

write  (46,4550)  ul,u2 
4550     format  (Ix,10x,f7._3,23x,f7.  3) 

write  (46,4600) 
4600     format (lx,/,lx, 'Side  Pier  Proportional  Limit ', 3x, 'Keel  Pier* 
+,'  Proportional  Limit') 

write (46, 4650)  plside,plside 
4650     format (lx,10x, f7.3,lx, 'kips/in2 ' 15x, f7 . 3, lx, 'klps/in2') 

write (46,4700) 
4700     format (lx,/,lx, 'Side  Pier  Vertical  Stiffness ', 3x, ' Side  Pier', 
+'  Horizontal  Stiffness') 

write (46,4750)  kvs,khs 
4750     format(lx,3x,fll.l,lx, 'kips/in' , llx, f 11 . 1, lx, 'kips/in') 

write  (46, 4775) 
4775    format (lx./,lx, 'Keel  Pier  Vertical  Stif fness ' , 3x, 
+'Keel  Pier  Horizontal  Stiffness') 

write (46,4780)  kvk,khk 
4780    format (lx,3x, fll. l,lx, 'kips/in' , llx, fll . 1, lx, 'kips/in') 

write  (46,4800) 
4800     format (lx,//, 20x, ' *  System  Parameters  and  Inputs  *') 

write (46,4850) 
4850     format (lx,/, lx, ' Input  Forcing  Function  is  Horizontal  Component', 
+'  of  the  1946  El  Centro') 

write  (46,4875) 
4875    format (lx, 20x, '  Earthquake  Acceleration  Time  History.') 

write (46,4995) 
4995     format (lx,/,lx, 'Vertical/Horizontal  Ground  Acceleration  Ratio' 
+,3x,'Data  Time  Increment') 

write (46, 4990)  amp,dtau 
4990    format(lx,10x,f6.3,t55,f6.3,lX, 'sec') 

write  (46,4900) 
4900    format (lx,/, lx, 'Gravitational  Constant ', 3x, 'Percent  System 
+  Damping') 

write  (46,4950)  gravity, zeta*100 . 
4950    format(lx,7x, f6.2,lx, ' in/sec2 ' , lOx, f6.2,lx, '%') 

write (46, 5000) 
5000    format (lx,/,25x, 'Mass  Matrix',/) 

do  5100  i=l,3 

write  (46,5050)  m(i,l)  ,m(i,2)  ,m(i,3) 
5050     format (lx, fl5.4,5x, fl5.4,5x, fl5.4) 
5100    continue 

write (46,5200) 
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5200     format (lx,/,25x, 'Damping  Matrix',/) 

do  5300  1=1,3 

write (46, 5250)  cx(i,l) ,cx(i,2) ,cx(i,3) 
5250     format (lx, fl5.4,5x, fl5.4,5x, fl5.4) 
5300    continue 

write (46, 5400) 
5400    format (lx,/,25x. 'Stiffness  Matrix',/) 

do  5500  1=1,3 

write  (46,5450)  k(i,l)  ,k(i,2)  ,k(i,3) 
5450     format ( lx , f 15 . 4 , 5x, f 15 . 4, 5x, f 15 . 4) 
5500    continue 

write (46, 5700) 
5700    format (lx,//) 

WRITE (46,6000) 

6000  FORMAT (IX, 'Undamped  Natural  Frequencies ' ,t35, 'Mode  #l',t50, 
+'Mode  #2',t65, 'Mode  #3') 

write (46,6001)  wl,w2,w3 

6001  format(lx,t31,f7.3,lx, 'rad/sec' , t46, f7 . 3, lx, 'rad/sec' ,t62,f7.3, 
+'  rad/sec') 

WRITE (46,6002) 

6002  FORMAT (IX, 'Damped  Natural  Frequencies ' ,t35, 'Mode  #l',t50, 
+'Mode  #2' ,t65, 'Mode  #3') 

WRITE (46,6500)  wl*sqrt (1-zeta* *2) , w2*sqrt (l-zeta**2) ,w3*sqrt (l-zeta**2) 
6500    format (lx,t 31, f7.3,lx, 'rad/sec' ,t46, f7.3,lx, 'rad/sec' ,t62,f7.3, 
+ '  rad/sec ' ) 
endlf 

write  (46,10500)  ampacc*100 
10500    format (lx,///, lx, 'For  Earthquake  Acceleration  of  ',f6.2,'  %  ' 
+,'of  the  El  Centro',/) 

write (46, 25000) 

25000  format (lx, 'Maxlmums/Failures * ,t26, 'X  (ins) ' , t36, ' Y  (ins)',t51, 
+'Theta  (rads) ' ,t65, 'Time  (sec)') 

write  (46, 25001) 

25001  format  (lx,  * '  ,t25,  ' '  ,t35,  ' ',t50, 

+  < •  /t64/  • •) 

write  (46,310)  maxx,timex 

310  format  (lx, '  Maximum  X' , t25, f9 .6, t65, f5. 2) 
write  (46,311)  maxy,timey 

311  format  (lx, '  Maximum  Y' , t35, f9.6,t65, f5. 2) 
write  (46,312)  maxt, timet 

312  format  (lx, '  Maximum  Rotation' ,t50, f9 .6, t65, f5. 2) 

if  (flagl.eq.l)  then 
flagl0=flagl0+l 

write  (46,313)  xl,yl,tl, timel 

313  format  (lx, ' Side  block  sliding'  , t25, f9.6, t35, f9.6, t50, f9.6, 
+t65,f5.2) 

endif 

if  (flag2.eq.l)  then 
flagl0=flagl0+l 

write  (46,314)  x2,y2, t2,time2 

314  format  (lx, 'Keel  block  sliding'  , t25, f9 .6, t35, f9.6, t50, f9 .6, 
+t65,f5.2) 

endif 

if  (flag3.eq.l)  then 


-136- 


flaglO=flaglO+l 

write  (46,315)  x3,y3, t3, time3 

315  format  (lx, 'Side  block  overturning'  , t25, f9.6, t35, f9.6, t50, f9.6, 
+t65,f5.2) 

end  if 

if  (flag4.eq.l)  then 
flaglO=flaglO+l 

write  (46,316)  x4,y4, t4.time4 

316  format  (lx, 'Keel  block  overturning'  , t25, f9.6, t35, f9.6, t50, f9.6, 
+t65,f5.2) 

endif 

if  (flagS. eq.l)  then 
flaglO=flaglO+l 

write  (46,317)  x5,y5, t5, tlme5 

317  format  (lx, 'Side  block  liftoff  , t25, f9.6,t35, f9.6, t50, f9 .6, 
+t65,f5.2) 

endif 

if  (flag6.eq.l)  then 
flaglO=flaglO+l 

write  (46,318)  x6,y6, t6, time6 

318  format  (lx, 'Keel  block  liftoff  , t25, f9.6, t35, f9.6, t50, f9.6, 
+t65,f5.2) 

endif 

if  (flag7.eq.l)  then 
flaglO=flaglO+l 

write  (46,319)  x7,y7, t7, time7 

319  format  (lx, 'Side  block  crushing*  , t25, f9.6,t35, f9.6, t50, f9.6, 
+t65,f5.2) 

endif 

if  (flag8.eq.l)  then 
flaglO=flaglO+l 

write  (46,320)  x8,y8, t8, time8 

320  format  (lx, 'Keel  block  crushing'  , t25, f9 .6,t35, f9 .6, t50, f9 .6, 
+t65,f5.2) 

endif 

if (flaglO.le.O)  then 
write (46, 11000) 
11000    format (lx,/,lx, 'No  failures  occurred.') 

if (counter .eq. 1 .0  .and.  flaglO.le.O)  then 

go  to  60000 

endif 

if (counter .eq.O .0)  then 

ampacmax=ampacc 

ampacc=ampacc+ . 1 

counter=l .0 

write  (6,*)  'I  am  in  the  secondary  looping  stage.' 

endif 

endif 

if  (ampacc . le .ampacmax)  go  to  20000 

if  (counter .eq.l .0)  then 

ampacc=ampacc- . 0 1 

else  if (counter .eq.O .0)  then 

ampacc=ampacc- . 1 

endif 
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go  to  10000 
20000   continue 
stop 
end 

subroutine  mult  (a,b,cyy) 

real  b(6,6) 

complex  a (6,6) ,cyy (6,6) ,d (6) 

do  100,1=1,6 

do  200,3=1,6 

do  300,k=l,6 

d(k)=a(i,k)*b(k,  j) 
300     continue 

cyy  (1 -  3)  =d  (1)  +d  (2)  +d  (3)  +d  (4)  +d  (5)  +d  (6) 
200     continue 
100     continue 

return 

end 

subroutine  multc  (a,b,cyx) 

complex  a(6,6)  ,b(6,6)  ,cyx(6,6)  ,d(6) 

do  100,  1=1,6 

do  200,  3=1,6 

do  300,  k=l,6 

d(k)=a(i,k)  *b(k,  j) 
300     continue 

cyx(i,  j)=d(l)+d(2)+d(3)+d(4)+d(5)+d(6) 
200     continue 
100     continue 

return 

end 
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APPENDIX  2 


C        NON-LINEAR  THREE  DEGREE  OF  FREEDOM  SYSTEM  RESPONSE 
C        USING  FOURTH  ORDER  RUNCE-KUTTA  METHOD 

integer  NN, 1 , mm, n, hull ,nsys, f laglO 

Integer  flagl, flag2, f lag3, flag4, flag5, flag6, f lag7, flag8 

real  beta, weight, h, Ik, gravity 

real  m(3,3) ,cx(3,3) ,k(3,3) ,sidearea,keelarea,plside 

real  bases ide,basekeel ,htside,htkeel,critl,crit2,crit3 

real  ac(6000) ,dtau,maxx,maxt,maxy,timex, timet 

real  rfl,rf2,rf3,hfl,hf2,hf3,  ampaccmass,  ampacmax 

rea 1  kvs , kvk , khs , khk , base , ht , counter 

real  timel , time 2 , time 3, time4, time5, time6, time7, time8 

real  x(6000)  ,t(6000)  ,y(6000)  ,XSCL(6) 

real  bbb,ccc,wl2 , wl , w2  2 ,w2,w32,w3, model, mode 3, capwidth 

real  mmxl,mmangl ,mmx3,mmang3,br , amp,plkeel ,crit4,ul,u2 

real  timey,mmininml,inmmmm2,mmmmm3,mmmmm4 

real  R,S,TAU,A(5)  ,B(5)  ,C(5)  ,D(5)  ,E(5)  ,F(5)  ,G(5)  ,HH(5) 

CHARACTER* 40  XLABEL, YLABEL, YYLABEL, YYYLABEL,DEC 


C  READ  IN  VESSEL  AND  DRYDOCK  DATA;  VESSEL  WEIGHT, KG, I (ABOUT  KEEL) 

C  TIME  INCREMENT  OF  DATA  POINTS, VERTICAL  STIFFNESS  OF  SIDE  AND 

C  KEEL  PIERS, HORIZONTAL  STIFFNESS  OF  SIDE  AND  KEEL  PIERS, 

C  GAVITATIONAL  CONSTANT, SIDE  BLOCK  BASE  AND  HEIGHT, 

C  KEEL  BLOCK  BASE  AND  HEIGHT, 

C  BLOCK-BLOCK  AND  BLOCK-HULL  FRICTION  COEFFICIENTS, 

C  SIDE  AND  KEEL  BLOCK'S  PROPORTIONAL  LIMIT, 

C  SIDE  PIER-VESSEL  CONTACT  AREA, KEEL  PIER-VESSEL  CONTACT  AREA, 

C  CAP  BLOCK  INCLINATION  ANGLE. 

read (44, *)  welght,h, Ik, dtau, kvs, kvk, khs, khk, gravity 

read(44,*)  baseside,  basekeel ,htside,htkeel ,ul,u2 

read (44, *)  br , amp,plside,plkeel , sidearea, keel area, zeta 

read (44,*)  hull,nsys 

write  (6,*)  'do  you  want  response  plots?  (y  or  n) ' 

read  (5, 15)  dec 
15      format  (a) 

do  12,1=1,3 

do  13,J=1,3 

m(l,j)=0.0 

k(i,j)=0.0 

cx(i,  j)=0.0 
13      continue 
12      continue 


CALCULATE  SYSTEM  PARAMETERS 
mass=weight/gravlty 
beta=asin(sqrt(br**2/(4*h**2) )) 
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m(l,  l)=mass 
111(1,  3)=h*mass 
m(2,2)=mass 
m (3, 1) =mass*h 
m(3,3)=Ik 

k(l,l)=(2*khs+khk) 
k(2,2)=(2*kvs+kvk) 
k(3,3)  =  ((0.5*kvs* (br**2 .0) ) - (weight *h) ) 

DETERMINE  NATURAL  FREQUENCIES  OF  SYSTEM 

bbb=-(m(l,l)  *k(3,3)+m(3,3)  *k  (1, 1) -m  (1,  3)  *k  (3, 1) -m  (3, 1)  *k  (1,  3) ) 
/(m(l,l)*m(3,3)-m(l,3)*m(3,l)) 

ccc=(k(l,l)*k(3,3)-k(l,3)*k(3,l))/(m(l,l)*m(3,3)-m(l,3)*m(3,l)) 


NATURAL  FREQ.  MODE  #1 

wl2= (-bbb-sqrt (bbb*  *  2 . 0-4*ccc) ) /2 
wl=sqrt (wl2) 

NATURAL  FREQ.  MODE  #2 

v22=k(2,2)/m(2,2) 
w2=sqrt (w22) 

NATURAL  FREQ.  MODE  #3 

w32= (-bbb+sqrt (bbb*  *2 . 0-4*ccc) ) /2 
w3=sqrt (w32) 

MODE  SHAPE  #1  &  #3 

model=(m(l,3)  *wl2-k(l,3)  )/(-m(l,l)  *vl2+k(l,l)) 
mode3=(m(l,3)  *w32-k (1, 3) ) / (-m (1, 1)  *v32+k(l,l)) 

DETERMINE  011,013,031,033 
mmxl=m (1, 1)  +m  (1 ,  3) /model 
mmangl=model*m (3, 1) +m (3, 3) 
mmx3=m (1 , 1)  +m  (1 ,  3) /mode3 
mmang3=mode3*m (3, 1) +m (3, 3) 
mmmmml=2*zeta*mmxl*wl 
mmmmm2=2  *  zeta  *  mmx3  *  w3 
mmmmm3=2 *zeta*mmangl * wl 
mmmmm4=2*zeta*mmang3*v3 


ex  (1,  3)=(mmmmml-mmmrom2)/(l/model-l/mode3) 
ex  (1, 1)  =mmmniml-  (ex  (1,  3)  /model) 
cx(2,2)=2*zeta*m(2,2) *w2 
ex (3, 1) = (mmmmm3-mmmmm4) / (model-mode3) 
ex (3, 3) =mmmmm3- (ex (3,1) *model) 


READ  IN  ACCELERATION  DATA 
do  300, n=l, 500 1,5 
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read  (45,*)  ac  (n) ,ac  (n+1)  , ac (n+2) , ac (n+3) ,ac (n+4) 
300     continue 
C       ESTABLISH  FAILURE  CRITERIA  AND  FLAGS 

critl=  mln(ul, (u2*cos (beta) +sin (beta) )/ 
+   (cos (beta) -u2*sin (beta) ) ) 
crit2=min  (ul,u2) 

crlt3=  (0 .66*baseside-12 .0) /htside 
crit4=basekeel/(6*htkeel) 
ampacc=l .0 
counter=0 .0 
ampacmax=0 . 0 
10000    continue 
flagl=0 
flag2=0 
flag3=0 
flag4=0 
flag5=0 
flag6=0 
flag7=0 
flag8=0 
flagl0=0 
maxx=0 . 0 
maxt=0.0 
maxy=0 . 0 
mm=0 

x(l)=0.0 
y(l)=0.0 
t(l)=0.0 
R=0.0 
S=0.0 
TAU=0.0 


C       IMPLEMENTATION  OF  EQUATIONS  OF  MOTION  INTO  THE 
C        RUNGE-KUTTA  FORMULUS 

do  301,1=1,6000 
DO  3000,11=0,4 

A(11)=0.0 

B(11)=0.0 

C(11)=0.0 

D(11)=0.0 

E(11)=0.0 

F(11)=0.0 

G(11)=0.0 

HH(11)=0.0 
3000     CONTINUE 
nun— uun+1 
DO  302,  NN=1,4 
IF(NN.EQ.l)  THEN 

FF=0.0 
ELSE  IF  (NN.EQ.2  .OR.  NN.EQ.3)  THEN 

FF=0.5 
ELSE  IF  (NN.EQ.4)  THEN 

FF=1.0 
ENDIF 
A(NN)=dtau*  (R+FF*D(NN-1) ) 


B (NN) =dtau* (S+FF*E (NN-1) ) 


-141- 


C (NN) =dtau* (TAU+FF*F (NN-1) ) 

D (NN) =dtau* ( (-ex (2 , 2) /m (2 , 2) ) * (R+FF*D (NN-1) ) - (k (2 , 2) /m (2 , 2) ) 
+  *(y(mm)+FF*A(NN-l))-amp*ampacc*ac(l)/2.54) 


G (NN) =dtau* ( (-ex (1 , 1) /m (1 , 1) ) * (S+FF*E (NN-1) ) - (ex (1 , 3) /m (1 , 1) ) 
+* (TAU+FF*F (NN-1) ) - (k (1 , 1) /m (1 , 1) ) * (x (mm) +FF*B (NN-1) ) -ampacc*ac (1 


)/ 


E(NN)  =  (m(l/l)*m(3,3)*G(NN)-m(l,3)*m(3/3)*HH(NN))/ 
+  (m(3,3)*m(l,l)-m(l,3)*m(3,l)) 

F  (NN)  =  (HH(NN)  -  (m  (3, 1) /m (3, 3) ) *E  (NN) ) 
302     continue 


C       DETERMINING  SYSTEM  RESPONSE 

y(mm+l)=y(mm)  +  (A(l)+2*A(2)+2*A(3)+A(4))/6 
x  (nnn+1)  =x  (mm)  +  (B  (1)  +2*B  (2)  +2*B  (3)  +B  (4)  )  /6 
t(mm+l)=t(mm)  +  (C(l)+2*C(2)+2*C(3)+C(4))/6 
R=R+(D(l)+2*D(2)+2*D(3)+D(4))/6 
S=S+(E(l)+2*E(2)+2*E(3)+E(4))/6 
TAU=TAU+ (F (1) +2*F (2) +2*F (3) +F (4) ) /6 

C       MAXIMUM  VALUES  FOR  TRANSLATIONS  AND  ROTATION 

If  (abs (x (mm) ) . gt .abs (maxx) )  then 

timex=dtau* (1-1) 

maxx=x (mml, 
endlf 
if  (abs (t (mm) ) .gt .abs (maxt) )  then 

tlmet=dtau* (1-1) 

maxt=t (mm) 
endif 
If  (abs (y (mm) ) .gt .abs (maxy) )  then 

tlmey=dtau* (1-1) 

maxy=y (mm) 
endlf 

C       CALCULATE  VERTICAL  AND  HORIZONTAL  FORCES  CAUSED  BY  VESSEL, 
C       TEST  FOR  FAILURE 

C        CALCULATE  FORCES  ON  SIDE/KEEL  BLOCKS 

rfl=kvs*  ((veight/k(2,2))-y(mm)-(br/2)  *t(mm)  ) 
r f 2=kvs* ( (weight A (2,2) )-y (mm) + (br/2) *t (mm) ) 
rf3=kvk*  ((welghtA(2,2)  )  -y  (mm)  ) 
hfl=khs* (x(mm)) 
hf2=khs* (x(mm)) 
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hf3=khk*x(mm) 
C       TEST  FOR  SIDE  BLOCK  SLIDING 

if  (flagl.eq.l)  then 

go  to  400 
else  if  (hfl.lt. 0.0. and. rfl.gt. 0.0. and. abs(hfl/rfl) .gt.critl)  then 

timel=  dtau* (1-1) 

flagl=l 
else  if  (hf2. gt. 0.0. and. rf2.gt. 0.0. and. abs(hf2/rf2) .gt.critl)  then 

timel=dtau* (1-1) 

flagl=l 
endlf 
xl=x (mm) 
yl=y (mm) 
tl=t (mm) 
400     continue 

C       TEST  FOR  KEEL  BLOCK  SLIDING 

if  (flag2.eq.l)  then 

go  to  410 
else  if  (rf3.gt.0.0.and.abs(hf3/rf3) .gt.crit2)  then 

time2=dtau* (1-1) 

flag2=l 
endif 
x2=x  (mm) 
y2=y  (mm) 
t2=t(mm) 
410  continue 
C       TEST  FOR  SIDE  BLOCK  OVERTURNING 

if  (flag3.eq.l)  then 

go  to  420 
else  if  (hfl. It. 0.0. and. rfl.gt. 0.0. and. abs (hfl/rfl) .gt.crit3)  then 

time3=  dtau*  (1-1) 

flag3=l 
else  if  (hf2.gt.0.0.and.rf2.gt.0.0.and.abs(hf2/rf2) .gt.crit3)  then 

time3=dtau*  (1-1) 

flag3=l 
endif 
x3=x (mm) 
y3=y (mm) 
t3=t (mm) 
420     continue 

C       TEST  FOR  KEEL  BLOCK  OVERTURNING 

if  (flag4.eq.l)  then 

go  to  430 
else  if  (rf3.gt.0.0.and.abs(hf3/rf3) .gt.crit4)  then 

time4=dtau* (1-1) 

flag4=l 
endif 
x4=x (mm) 
y4=y (mm) 
t4=t (mm) 
430     continue 

C       TEST  FOR  SIDE  BLOCK  LIFTOFE 
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If  (flag5.eq.l)  then 

go  to  440 
else  If  (rfl.lt. 0.0  .or.  rf2.lt. 0.0)  then 

tlme5=dtau* (1-1) 

flag5=l 
endlf 
x5=x (mm) 
y5=y  (mm) 
t5=t (mm) 
440     continue 

C       TEST  FOR  KEEL  BLOCK  LIFTOFF 

if  (flag6.eq.l)  then 

go  to  450 
else  if  (rf3.lt. 0.0)  then 

tlme6=dtau* (1-1) 

flag6=l 
end  if 
x6=x (mm) 
y6=y (mm) 
t6=t (mm) 
450     continue 

C       TEST  FOR  SIDE  BLOCK  CRUSHING 

if  (flag7.eq.l)  then 

go  to  460 
else  if  (rfl.gt.0.0  .and.  (r f 1/sldearea) .gt .plside)  then 

flag7=l 

time7=dtau* (1-1) 

else  if  (rf2.gt.0.0  .and.  (r f2/sidearea) .gt. plside)  then 
flag7=l 

time7=dtau* (1-1) 
endif 
x7=x (mm) 
y7=y (mm) 
t7=t (mm) 
460     continue 

C       TEST  FOR  KEEL  BLOCK  CRUSHING 

if  (flag8.eq.l)  then 

go  to  470 
else  if  (rf3.gt.0.0  .and.  (r f3/keelarea) .gt .plkeel)  then 

flag8=l 

time8=dtau* (1-1) 
endlf 
x8=x (mm) 
y8=y (mm) 
t8=t (mm) 
470     continue 

301     continue 


C       PLOT  RESULTS 

go  to  999 
60000   continue 
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998 
999 


if (DEC.EQ. 'N')  THEN 

write (6,*)  'I  am  finishing.' 

GO  TO  998 
endif 

vrite(6,*)  'I  am  plotting. ' 
XSCL(1)=0.0 
XSCL (2) =30.0 
XLABEL='TIME  IN  SECONDS' 
YLABEL=' ROTATION  IN  RADIANS' 

YYLABEL= 'RELATIVE  HORIZONTIAL  DISPLACEMENT  IN  INCHES' 
YYYLABEL= 'RELATIVE  VERTICAL  DISPLACEMENT  IN  INCHES' 
CALL  QPICTR  (X, 1, 3000,QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) , 
QYLAB(YYLABEL) , QLABEL (4)) 

CALL  QPICTR  (T, 1, 3000, QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) , 
QYLAB (YLABEL) , QLABEL (4) ) 

CALL  QPICTR  (Y, 1, 3000, QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) , 
QYLAB (YYYLABEL) , QLABEL (4)) 
go  to  20000 
CONTINUE 


if (ampacc .eq. 1 .0)  then 


4000 
4050 
4100 
4150 
4200 

4250 
4300 

4350 

4400 
4450 
4470 

4475 
4500 

4550 
4600 

4650 


Hull  ',I3,lx, '**') 
Ship  Parameters  *') 


Drydock  Parameter  * ' ) 


write  (46,4000)  nsys 

format (lx,/,28x, '****  System  ', 12 , lx, '****' ) 

write (46,4050)  hull 

format (lx,/,30x, '** 

write  (46,4100) 

format (lx,//,28x, '* 

write  (46,4150) 

format (lx,/,5x, 'Weight' ,8x, 'Moment  of  Inertia' ,9x, 'K.G. ' ) 

write  (46, 4200)  weight, Ik, h 

format (lx, f9.1,lx, 'kips' , lx, f 11 . 1 , lx, "kips-ln-sec2 ' , 
+3x,f6.1,lX, 'ins') 

write  (46,4250) 

format (lx,//, 26x, 

write  (46, 4300) 

format(lx,/,lx, 'Side  Block  Height ', 3x, ' Side  Block  Width', 
+3x,'Keel  Block  Height ', 3x, 'Keel  Block  Width') 

write (46,4350)  htslde,baseside,htkeel ,basekeel 

format (2x, f6 . 1 , lx, ' ins ' , llx, f6 . 1 , lx, ' ins ' , llx, f6 .1, lx, ' ins ' , 
+9x,f6.1,lx, 'ins') 

write (46, 4400) 

format (lx,/,lx, 'Side-to-Side  Pier  Distance') 

write (46,4450)  br 

format (lx,t7, f6.1,lx, 'ins') 

write (46,4470) 

format (lx,/, '  Total  Side  Pier  Contact  Area' 
+,3x, 'Total  Keel  Pier  Contact  Area') 

write (46,4475)  sldearea,keelarea 

format (lx,8x,f 11. l,lx, 'in2' , 14x, f 11 . 1, lx, 'in2') 

write (46, 4500) 

format (lx,/, lx, 'Block-on-Block  Eriction  Coef f * ,3x, 'Hull-on-Block 
+  Eriction  Coeff) 

write (46, 4550)  ul,u2 

format (lx,10x, f7 . 3, 23x, f7. 3) 

write  (46,4600) 

format (lx,/,lx, 'Side  Pier  Proportional  Limit '. 3x, 'Keel  Pier' 
+,'  Proportional  Limit') 

write  (46, 4650)  plside,plside 

format (lx,10x, f7.3,lx, 'kips/in2 ' 15x, f7. 3, lx, ,kips/ln2') 
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write (46,4700) 
4700    format (lx,/,lx, 'Side  Pier  Vertical  Stiffness ', 3x, 'Side  Pier', 
+'  Horizontal  Stiffness') 

write  (46,4750)  kvs,khs 
4750     format(lx,3x,fll.l,lx, 'kips/in' , llx, fll . 1, lx, 'kips/in') 

write (46,4775) 
4775     format (lx,/,lx, 'Keel  Pier  Vertical  Stlf fness ' , 3x, 
+'Keel  Pier  Horizontal  Stiffness') 

write (46,4780)  kvk,khk 
4780     format(lx,3x,fll.l,lx, 'kips/in' , llx, fll . 1, lx, 'kips/in') 

write (46, 4800) 
4800    format (lx,//, 20x, *  *  System  Parameters  and  Inputs  **) 

write (46,4850) 
4850    format (lx,/, lx, ' Input  Forcing  Function  is  Horizontal  Component', 
+'  of  the  1946  El  Centro') 

write (46, 4875) 
4875    format (lx, 20x, '  Earthquake  Acceleration  Time  History.') 

write  (46,4995) 
4995    format (lx,/,lx, 'Vertical/Horizontal  Ground  Acceleration  Ratio' 
+,3x,'Data  Time  Increment') 

write (46, 4990)  amp,dtau 
4990     format (lx,10x,f6.3.t55,f6.3,lX. 'sec') 

write  (46,4900) 
4900    format (lx,/, lx, 'Gravitational  Constant ', 3x, 'Percent  System 
+  Damping') 

write (46,4950)  gravity, zeta*100 . 
4950     format (lx,7x,f 6.2, lx, 'in/sec2' , lOx, f6 . 2, lx, '%') 

write  (46,5000) 
5000     format (lx,/,25x, 'Mass  Matrix',/) 

do  5100  1=1,3 

write  (46,5050)  m(i,l),m(i,2),m(i,3) 
5050    format (lx, fl5.4,5x, f 15.4, 5x, f 15.4) 
5100     continue 

write (46,5200) 
5200     format (lx,/,25x, 'Damping  Matrix',/) 

do  5300  1=1,3 

write(46,5250)  cx(i,l) ,cx(i,2) ,cx(i,3) 
5250     format (lx, fl5.4,5x, fl5.4,5x, fl5 .4) 
5300    continue 

write (46, 5400) 
5400    format (lx,/,25x, 'Stiffness  Matrix',/) 

do  5500  1=1,3 

write  (46, 5450)  k(i,l),k(i,2),k(i,3) 
5450     format (lx , f 15 . 4 , 5x , f 15 . 4, 5x, f 15 . 4) 
5500    continue 

write  (46,5700) 
5700    format (lx,//) 

WRITE (46,6000) 

6000  FORMAT (IX, 'Undamped  Natural  Frequencies ' ,t35, 'Mode  #l',t50, 
+'Mode  #2',t65, 'Mode  #3') 

write (46,6001)  wl,w2,w3 

6001  format (lx,t31, f7. 3, lx, 'rad/sec' , t46, f7 . 3, lx, 'rad/sec' ,t62, f7.3, 
+'  rad/sec') 

WRITE (46,6002) 

6002  FORMAT (IX, 'Damped  Natural  Frequencies  * ,t35, 'Mode  #l',t50, 
+'Mode  #2*,t65, 'Mode  #3') 

WRITE (46,6500)  wl*sqrt (l-zeta**2) ,w2*sqrt (l-zeta**2) ,w3*sqrt (l-zeta**2) 
6500     format(lx,t31,f7.3,lx, 'rad/sec' , t46, f7.3, lx, 'rad/sec' ,t62,f7.3, 
+'  rad/sec') 
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endif 

write  (46, 10500)  ampacc*100 
10500    format (lx,///,lx, 'For  Earthquake  Acceleration  of  ',f6.2,'  %  * 
+  ,*of  the  El  Centro',/) 

write  (46, 25000) 

25000  format (lx, 'Maximums /Failures ' ,t26, 'X  (ins) ' ,t36, 'Y  (ins)',t51, 
+'Theta  (r ads) ' ,t65, 'Time  (sec)') 

write(46, 25001) 

25001  format  (lx,  ' '  ,t25,  ' '  ,t35,  ' '  ,t50, 

+  • '  ,  t64,  ' ' ) 

write  (46,310)  maxx,timex 

310  format  (lx, *  Maximum  X" , t25, f9.6, t65, f5. 2) 
write  (46,311)  maxy,tlmey 

311  format  (lx, '  Maximum  Y' , t35, f9 .6,t65, f5. 2) 
write  (46,312)  maxt, timet 

312  format  (lx, '  Maximum  Rotation' ,t50, f 9 .6, t65, f5. 2) 

if  (flagl.eq.l)  then 
flagl0=flagl0+l 

write  (46,313)  xl.yl, tl,timel 

313  format  (lx, ' Side  block  sliding'  , t25, f9 .6, t35, f9.6, t50, f9 .6, 
+t65,f5.2) 

endif 

if  (flag2.eq.l)  then 
flagl0=flagl0+l 

write  (46,314)  x2,y2, t2, time2 

314  format  (lx,'Keel  block  sliding'  , t25, f9.6,t35, f9.6, t50, f9 .6, 
+t65,f5.2) 

endif 

if  (flag3.eq.l)  then 
flagl0=flagl0+l 

write  (46,315)  x3,y3, t3,time3 

315  format  (lx, ' Side  block  overturning'  , t25, f9.6,t35, f9.6,t50, f9.6, 
+t65,f5.2) 

endif 

if  (flag4.eq.l)  then 
flagl0=flagl0+l 

write  (46,316)  x4,y4, t4, time4 

316  format  (lx, 'Keel  block  overturning'  , t25, f9.6,t35, f9.6, t50, f9 .6, 
+t65,f5.2) 

endif 

if    (flag5.eq.l)    then 
flagl0=flagl0+l 

write    (46,317)    x5,y5, t5,tlme5 

317  format    (lx, 'Side  block    liftoff    , t25, f9.6, t35, f9.6, t50, f9.6, 
+t65, f5.2) 

endif 

if  (flag6.eq.l)  then 
flagl0=flagl0+l 

write  (46,318)  x6,y6,t6,time6 

318  format  (lx, 'Keel  block  liftoff  ,t25, f9 .6, t35, f9.6, t50, f9 .6, 
+t65,f5.2) 
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endif 


319 


if  (flag7.eq.l)  then 
flaglO=flaglO+l 

write  (46,319)  x7,y7, t7, time7 
format  (lx, 'Side  block  crushing'  ,t25, f9 .6, t35, f9.6, t50, f9.6, 
+t65,f5.2) 
endif 


320 


if  (flag8.eq.l)  then 
flagl0=flagl0+l 

write  (46,320)  x8,y8, t8, time8 
format  (lx, 'Keel  block  crushing'  , t25, f9 .6, t35, f9.6, t50, f9.6, 
+t65,f5.2) 
endif 


if (flaglO.le.O)  then 

write (46, 11000) 
11000    format (lx,/,lx, 'No  failures  occurred.') 

if (counter .eq.l .0  .and.  flaglO.le.O)  then 

go  to  60000 

endif 

if (counter .eq.O .0)  then 

ampacmax=ampacc 

ampacc=ampacc+ . 1 

counter=l .0 

write  (6,*)  'I  am  in  the  secondary  looping  stage. 

endif 

endif 

if (ampacc . le . ampacmax)  go  to  20000 

if  (counter .eq. 1 .0)  then 

ampacc=ampacc- .01 

else  if  (counter .eq.O .0)  then 

ampacc=ampacc- . 1 

endif 

go  to  10000 
20000   continue 

stop 

end 
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APPENDIX  3 
MODAL  ANALYSIS  OF  THE  TWO  AND  THREE  DEGREE  OF  FREEDOM  SYSTEMS 

Since  the  vertical  equation  of  motion  is  uncoupled  from 
the  other  equations,  the  three  degree  of  freedom  system  need 
only  to  be  analyzed  to  obtain  damping  coefficients  for  both 
degree  of  freedom  mathematical  models.   Also,  the  maximum 
response  of  the  two  systems  will  be  determined  using  response 
spectrum  analysis  with  participation  factors  described  in 
Section  6.4. 

The  three  degree  of  freedom  equations  of  motions, 
undamped,  as  shown  in  Equations  4.4.11a,  b  and  c  are 

mn  x  +  m^3  6  +  k]^  x  =  -m11  x„ 

m22  Y  +  k22  y  =  ~m22  Y"g 
and 

m33e     +  m31   x  +   k33e     =   -m31   xg 

where  m11  =  m22  =  M 

m13  =  m31  =  MKG 

m33  "  xk 

kll  =  2Ksh  +  Kkh 

k22  =  2Ksv  +  Kkv 

k33  =  (B2/2)  Ksv  -  W  KG. 
To  perform  modal  analysis,  consider  the  free  vibration  system 

mn  x  +  m13  6  +  kn  x  =  0 

m22  y  +  k22  y  =  ° 

and 

m33  6+m31x+k33  6=0. 
Assume  the  system  response  is  in  the  form 
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x  =  4>i   sin(wt  +  P  ) 

y  =  <J>2  sin(cot  +  p ) 
and  G  =  $3  sin  (cot  +p )  . 
Now,    the  equations  of  motion  are 

mll(""<t>iw2   sin  (cot  +  p  ))  +  m13    (-tf)  3aj2sin(w  t  +  p)) 

+  kn    (<J>  i   sin  (cot  +  p  ))  =  0 


and 


m22  (-4)2w2  sin  (cot  +  P))  +  k22  (^2  sin  (cot  +  P))  =  0 


m33(~<}>3  co2  sin  ^t  +  p  ))  +  m31  (-42  co2  sin  (cot  +p)) 

+  k33  (<j>3  sin  (cot  +  p  ))  =  0. 
The  trivial  solution  to  the  above  set  of  equations  is 
sin  (cot  +  p)  =0.   Assuming  sin  (cot  +p)  =  0,  then 


kirmnw' 


-m31w 


0  -m,  -.    co 

k22-m22a32  0 

0     k33-m33W  . 


r*r 

"o" 

*2 

= 

0 

_*3J 

_0_ 

(A. 3.1) 


These  equations  are  solved  for  w  by  setting  the  determinant  of 
the  first  matrix  equal  to  0.   However,  the  second  row  of  the 
matrix  is  uncoupled.   Thus  one  solution  for   is 

co   =  k22/in22- 
Now,  the  determinant  equation  reduces  to 


4     2 

CO      CO 


/mllkll+m33k33\  +    kllk33       =  Q 
\mllm33"m31m3l/  "  mllm33-m13m31 

Using  the  quadratic  formula,  the  other  values  of  co 2  are 

(A. 3. 2) 


-Y+UJ-, 


2        -bVb  -4c 

CO   =  
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where   b  = 


/mllkll+m33k33\ 
\mllm33"m13m3l/ 


and 


11  33 
c  = 


mllm33-ml3m31 


Using  parameters  found  in  Section  5  for  the  eleven  representa- 
tive vessel-drydock  systems,  the  natural  frequencies w ±,   u2, 
and  003  for  each  system  is  calculated  and  listed  in  Table  A3 . 1 

The  mode  shapes  are  determined  with  the  relationship 

(kll   "  mll  w2)    61   "  ml^2  ^3   =   0 

°r        *1  m13    w 

*3         kll"mllco2       '  1      3 

With  mode  shapes  determined,  the  equation  or  motion  can  be 

rewritten,  including  coupled  damping 

m±i   x  +  m13  9  +  c^  x  +  c13  9  +  kn  x  =  0 

m22  y  +  c22  y  +  k22  y  =  ° 

and 

m33  0+  m31  x  +  c33  6  +  c31  x  +  k33   =  0. 
At  natural  frequencies  wlf  w2  andw3  the  coupled  three  degree 
of  freedom  system  acts  as  a  single  degree  of  freedom  system 
with  the  following  relationship: 
(<J>l/<J>3)  x  =  e 

(<f>l/<i>3)    *  =  § 
and      C^i/^3)    x  =  9  . 

Using  these  relationships,  the  three  degree  of  freedom 
equations  of  motion  at  the  three  system  natural  frequencies 
are 
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[Ill    +    (<J>3/4>l)»13]    x   +    tcH    +    W3/h)c3l]    x   +    kll    x   =    ° 

•  ■  • 

m22  y  +  c22  y  +  k22  y  =  ° 

and 

•  •  • 

[(4>l/4>3)m31   +  m33]e    +    [(*!/<()  3 )c31   +   c33]9    +   k339     =    0. 

The  damping  coefficients  for  the  problem  are  5%  of  the 

critical  damping  for  each  mode, 

[Cn  +  (<t>3/^i)c13]  =  2  CMX  a)n  ,  con  =  1,3       (A3. 3a) 
c22  =  2£  m22  0)2  (A3. 3b) 

and 

[(4>l/<J>3)h  c31  +  C33]  =   2  M3wn/   (Jn=  1,3      (A3. 3c) 
where  ►   =  percentage  of  critical  damping,  0.05 

Ml  =  mll  +  ($3/$  l)n  m13 
and  M3  =  m33  +  (4>  i/(j)  3 )  n  m3 1  • 

Substituting  known  values  and  rearranging  Equations  A3. 3a  and 
A3. 3c  yields  the  following  four  equations: 
Mode  1,        Cn  +  (<J>  5/^1)1  c13  =  2  C  Mi  w  ]_ 

(4>l/<J>  3)i  c31  +  c33  =  2  £M3  a)  ! 
Mode  2,        Cn  +  (^/(J)!^  c13  =  2  5  Mx  w  3 

(*  l/*3)3  c31  +  c33  =  2  ?M3  w3 
or  clx  +  ($  3/t>  i)  !  c13  =  2  £  MX  w  1 

cn  +  (4>  3/4>  1)  3  c13  =  2  £  Kt   w  3 
and  (^i/^3)i  c31  +  c33  =  2?M3o)1 

(45  i/4>  3)3  c31  +  c33  =  2  £  M3  w  3 
Using  the  above  set  of  equations  and  equation  A3. 3b,  damping 
coefficients  c11(  c^3/  c22/  c3l/  anc*  c33  are  determined  and 
listed  in  Table  A3 . 2  for  the  eleven  vessel-drydock  configura- 
tions.  Note  that  c^3  equals  to  c3i«  In  the  two  degree  of 
freedom  case,  the  damping  coefficients  c^,  ci2'  c21'  anc*  c22 
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are  equal  to  those  in  the  three  degree  problem,  i.e. 
Two  Degree  Case  Three  Degree  Case 

cll  cll 

c12  c13 

c21  c31 

c22  c33 

Using  the  response  spectrum  method  and  modal 

participation  factors,  the  maximum  response  of  the  two  and 

three  degrees  of  freedom  systems  can  be  calculated.   Starting 

with  the  forced  vibration  equation 

!>H  +    (<J>3/*l)n  m13]    x  +  ci  x  +  kn  x  =  -mn  xg 

or         ••         cl*   +  kUx  "mll 

x   +  -p—  -.=   —  — .    xg 


^11+(^/^)nm13l       fm11+(*3/«1)m13i 


Equations  6.3.1,  6.3.2  and  6.3.3  predict  a  maximum  response 
for  a  5%  critically  damped  system  for  a  given  mode  as 

m11 

(x         )       =    11.62*   > — t- — T- — > =.      for   u),,oo-,    <    2.24 

max'n  fmH    f^/^)  nm13l 


or    (x        )      =    26,03    *    T J:2-    /x    s    -     -,  for    2.24    <    w,  ,u)-    <    12.74 

max   n  oo  |mn  +  ((l)V<l>i)    mi  ol  1      j 


KT^VhV^ 


m 

else    (x        )      =    331:53    *   T rri    n    ; =i     for   uj.  ,    w0    >    12.74 

max  n      2      fni.. ,  +  (c(u/<j),  )  m,  ri       1'   3 


CO 


3] 


Similarly  for  rotational  response,  the  equation  of  motion  is 
[m33  +  (cf,  i/4>3)n  m31]  e  +  c36+  k33  6  =  -m31  xg 


or 


..       c36+k33e  "m31  " 

[m33+  (<D1/*3)nm3  J  ~  [m33-f(^1/(J)3)nm31]   xg 
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and  the  maximum  response  for  given  mode  is 

m 
<9max»n  =  1X-62  [m33+(  ±/   ^m^J  for  V  "3  <  2'24 

or    (emax»n  "  ~  *  tn^+^/^^J  for  2"24  <  V  "3  <  12"74 

331.53        m31 
else  (6    )   =  = —  *  -for  cu   to   >  12.74 

^^      "2      ^'VVn^      '    3 


Since  being  uncoupled,  the  vertical  equation  of  motion  is 

m22  y  +  c2  ^  +  k22  y  "  ~m22  yg 

or    ••    c2y+k22y 

y  +  =  -V   =  -AMP  x 

*  m22        ->g  g 

where  AMP  =  ratio  of  vertical/horizontal  El  Centro  earthquake 

accelerations . 
The  maximum  vertical  response  is 

(ymax)2  =11.62   *  AMP     for  ^2  <  2/24 
or 

(ymax)2  =  26,03   *  AMP     for  2.24  <  ^2    <    12.74 

CO 

else 

(ymax)2  =  331_53  *  AMP     for  ^2  >  12.74. 

CO 

Now,  the  maximum  response  for  each  mode  of  the  x  and  e 
equations  and  for  the  uncoupled  y  equation  are  calculated. 

The  configuration  of  the  three  degree  of  freedom  system 
at  any  time  is  a  superposition  of  the  two  coupled  natural  mode 
shapes  along  with  the  independent  vertical  maximum  response. 
The  absolute  maximum  response  to  a  given  earthquake  is  the 
numerical  summation  of  the  maximum  response  of  each  mode  shape 
times  its  respective  participation  factor.   The  general 
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formula  for  each  modal  participation  factor  is 

r  n  =ZMr(J)  rn/  Mr  (J)  rn 
where     r n  =  modal  particpation  factor  for  the  nth  mode 
r~  refers  to  which  equation  of  motion  is  being 
considered.   In  the  three  degree  of  freeddom 
case,  r  =  1  or  3. 
Using  these  particpation  factors,  the  maximum  response  of  the 

three  degree  of  freedom  system  should  be  no  greater  than 

r  r 

xmax  =    (xmax)l      1  +    (xmax) 3      3' 

Ymax  =    (^max)2 

and     Qmax  =    (9max)l     fi+    (6max)3r     3- 

These  maximum  response  are  determined  and  listed  in  Table  A3 . 3 

for  the  eleven  vessel-drydock  configurations.   Once  again  due 

to  the  uncoupled  vertical  response,  the  maximum  response  of  x 

and  e  in  the  two  degree  of  freedom  is  identical  to  that  of  the 

three  degree  of  freedom  case. 
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SYSTEM 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 


TABLE  A3.1 

Three  Degree  of  Freedom  Vessel-Drydock  System 
Natural  Frequencies 
(rad/sec) 


wl 
5.454 

3.883 

4.013 

2.819 

2.848 

9.690 

7.953 

6.976 

5.863 

6.857 

6.177 


(0  2 
39.071 
36.823 
29.468 
27.773 
35.737 
45.450 
41.322 
39.247 
32.829 
34.832 
33.293 


0)  3 
37.230 
34.367 
23.791 
21.511 
32.742 
40.863 
36.351 
34.045 
29.060 
29.921 
28.010 


NOTE:   Two  Degree  of  Freedom  System  Natural  Frequencies 
correspond  to  w^  and  co3. 
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TABLE  A3. 2 

Three  Degree  of  Freedom  Vessel-Drydock  System 
Damping  Coefficients 
(kim-sec/inch) 


SYSTEM 

cll 

c22 

c33 

c13'  c31 

1 

65.398 

166.826 

1322882.125 

3920.416 

2 

58.612 

157.226 

938583.250 

2874.701 

3 

42.452 

125.824 

976657.438 

2829.964 

4 

37.308 

118.583 

682485.813 

2053.786 

5 

54.303 

152.588 

687504.500 

2159.187 

6 

188.657 

445.695 

7879946.000 

17130.227 

7 

167.100 

405.212 

6430934.500 

14269.135 

8 

155.786 

384.866 

5622059.000 

12660.506 

9 

43.963 

116.475 

1196244.250 

3340.509 

10 

33.833 

86.437 

836179.438 

2408.758 

11 

31.639 

82.618 

751062.813 

2185.568 

NOTE:   Two  Degree  of  Freedom  System  Damping  Coefficients 
correspond  to  c^,  C13/  c31  anc*  c33# 
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TABLE  A3. 3 

Three  Degree  of  Freedom  Vessel-Drydock  System 
Maximum  Response  Using  Response  Spectrum  Analysis 


SYSTEM 

fins. ) 

(rhs.) 

8 
(rads. ) 

1 

0.33275 

0.10859 

0.01646 

2 

0.27141 

0.12225 

0.02348 

3 

0.65142 

0.19089 

0.02350 

4 

0.53788 

0.21491 

0.02350 

5 

0.21468 

0.12980 

0.03202 

6 

0.45271 

0.08025 

0.00683 

7 

0.46833 

0.09708 

0.00866 

8 

0.46597 

0.10762 

0.01009 

9 

0.64366 

0.15381 

0.01525 

10 

0.65864 

0.13662 

0.01342 

11 

0.67467 

0.14955 

0.01522 

NOTE:   Two  Degree  of  Freedom  System  Maximum  Response 
correspond  to  x  and  6. 
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